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BALLISTIC  EVALUATION  OF  PROPELLANTS 
IN  MICRO- MOTORS 

ABSTRACT 

Development  of  a  10-gm  micro-motor  for  ballistic 
evaluation  of  novel  propellantB  was  initiated  by  this  Division  in  early 
I960.  Several  thousand  micro-motors  containing  a  variety  of  propellant 
formulations  were  fired  demonstrating  that  reliable  precise  ballistic 
data  could  be  obtained  if  proper  experimental  techniques  were  used. 

An  extensive  series  of  static  tests  showed  that  specific  • 
impulse  could  be  scaled  when  corrections  were  made  for  heat-loss  and  . 
two-phase-flow  effects.  These  scaling  factors  could  be  obtained  from  a 
series  of  micro-motor  tests  requiring  as  little  as  2  pounds  of  propellant. 

This  report  1b  a  summary  of  all  previous  reports  on  the 
development  of  micro-motors  and  small-scale  testing  and  it  is  directed 
toward  those  groups  which  might  wish  to  use  the  Rohm  &  Haas  micro¬ 
motor  system  for  propellant  evaluation  purposes.  Accordingly,  the 
areas  of  system  design  features,  motor  preparation  and  assembly,  ignition, 
data  acquisition,  data  reduction,  scaling  effects,  scaling-factor  determination, 
and  specifio-impulse  prediction  have  been  treated  in  considerable  detail. 
Sufficient  information  is  provided  to  allow  interested  agencies  to  duplicate 
the  hardware  system,  conduct  static  tests,  obtain  valid  ballistic  data  on 
scaling  factors,  and  make  specific-impulse  predictions. 
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BALLISTIC  EVALUATION  OF  PROPELLANTS 
IN  MICRO- MOTORS 


1 .  Introduction 

Over  the  past  several  years  there  has  been  a  general 
trend  in  propellant  research  toward  the  use  of  smaller  motors  for  ballistic- 
evaluation  purposes.  Many  factors  have  contributed  to  this,  including 
(a)  a  large  number  of  new  contractors  in  the  field  of  propellant  research, 
most  of  whom  have  an  upper  restriction  on  propellant  quantity,  (b)  increasing 
toxicity  and/or  sensitivity  of  new  formulations,  and  (c)  limited  availability 
or  excessive  costs  of  novel  propellant  ingredients.  These  factors  make 
ballistic  evaluation  of  propellants  in  small  motors  either  highly  desirable 
or  mandatory,  and  have  caused  many  groups  to  revaluate  their  ideas 
concerning  motor  si2e  and  ballistic  measurements. 

Early  in  I960  the  Redstone  Arsenal  Research  Division  of 
Rohm  &  Haas  Company  initiated  the  development  of  a  ten- gram  micro* 
motor  for  solid  propellant  static -test  purposes.  The  resulting  motor  was 
a  small-scale  model  (.75-in,  I.  D.  X  1.5-in.  long)  of  the  Division's  standard 
6C5-1X.4  test  motor.  Additional  ,75-in. -I.  D,  micro-motors  in  lengths 
of  2.5  and  3.5  inches  were  subsequently  developed.  Static  tests  of  several 
thousand  micro-motors  containing  a  variety  of  propellant  formulations 
demonstrated  that  reliable,  precise  ballistic  data  could  be  obtained  if 
proper  experimental  techniques  were  used. 

Concurrent  with  the  r'icro-motor  development,  a  parallel 
investigation  of  the  specific-impulse  scaling  phenomena  of  solid  propellants 
was  being  conducted  by  this  Division.  An  extensive  series  of  static  tests 
showed  that  specific  impulse  could  be  scaled  when  corrections  were  made 
for  heat-loss  and  two-phase-flow  effects.  Specifically,  it  was  found  that 
experimental  scaling  factors  relating  to  heat  flux  and  condensed-phase 
particle  diameter  could  be  extracted  from  data  obtained  by  static  testing 
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a  propellant  in  a  series  of  small-motor  configurations,  that  these  scaling 
factors  could  be  used  to  accurately  predict  the  specific  impulse  of  the 
propellant  in  any  other  untested  motor  configuration,  and  that  these 
scaling  factors  could  be  obtained  from  a  series  of  micro-motor  tests 
requiring  as  little  as  two  pounds  of  propellant. 

Most  recently,  the  micro-motor  hardware  and  casting 
equipment  were  redesigned  to  allow  completely  remote  operation,  and 
remote  casting  and  firing  facilities  were  designed  and  constructed.  These 
modifications  allow  evaluation  of  very  hazardous  materials. 

This  report  is  a  summary  of  all  previous  reports  on  the 
development  of  micro-motors  and  small-scale  testing  and  is  directed 
toward  those  groups  who  might  wish  to  use  the  Rohm  &  Haas  micro-motor 
system  for  propellant  evaluation  purposes.  Accordingly,  the  areas  of 
system  design  features,  motor  preparation  and  assembly,  ignition,  data 
acquisition,  data,  reduction,  scaling  effects,  scaling-factor  determination, 
and  specific  impulse  prediction  have  been  treated  in  considerable  detail. 
Sufficient  information  is  provided  to  allow  interested  agencies  to  duplicate 
the  hardware  system,  conduct  static  tests,  obtain  valid  ballistic  data 
and  scaling  factors,  and  make  specific-impulse  predictions.  Details  of 
previous  scaling  data  and  prediction  results  for  specific  propellants  may 
be  obtained  from  the  material  listed  in  the  Bibliography. 

?..  Design  Features 

2.1  General 

The  standard  Rohm  &  Haas  motor  designation  used  through¬ 
out  this  report  is  as  follows:  reading  from  left  to  right,  the  first  number 
represents  motor  I.  D.  ,  the  letter  group  designates  grain  design,  the  next 
number  represents  grain  I.  D.  ,  and  the  last  number  represents  motor 
chamber  length.  As  an  illustration  the  .75C.50-1.5  motor  has  a  case 
I.D.  of  0.75  in.  ,  a  cylindrical  grain  design,  a  grain  I.  D,  of  0.50  in.  ,  and 
a  motor  case  length  of  1.5  in. 
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2.2  Motors 

2.2.1  Thread-Type  .75C.5Q  Motor* 

The  threaded  .7  5C.50  motor  (Fig.  1)  is  attached  to  its 
associated  hardware  with  machine  threads  and  was  designed  for  non- remote 
assembly.  The  motor  bodieu  are  1.000-in.  O.  D.  X  0.750-in.  1.  D.  mild 
steel  tubes  with  1-1/8  -  12  NF  male  threads  on  each  end.  The  standard 
motor  lengths  are  1.5,  2.5,  and  3.5  in.  The  nozzle  is  mild  steel  with  a 
graphite  insert  and  is  threaded  on  the  exit  end  for  attachment  of  a  nozzle* 
closure  retainer.  The  other  end  of  the  motor  attaches  to  a  special  pressure 
transducer  and  squib-lead  sealing  assembly  which  is  in  turn  attached  to 
a  load  cell.  Squib  leads  are  sealed  by  a  gland  nut  and  Teflon  packing. 

Several  motors  may  be  connected  together  with  threaded 
couplings  to  obtain  greater  lengths.  Threaded  connections  provide  adequate 
sealing  for  all  joints  in  the  motor  system. 

Motors  nave  been  hydrostatically  tested  to  13,  000  psl 
without  failure, indicating  that  a  working  pressure  of  10,  000  psi  can  be 
considered  safe.  However,  an  all-steel  nozzle  is  required  for  pressures 
about  6000  psi  to  prevent  failure  of  the  graphite  insert. 

2.2.2  Clamp- Type  .75C.50  Motor 

The  clamp  type  .75C.50  motor  (Fig,  2)  was  designed  for 
remote  casting,  assembly,  and  firing  but  is  also  used  for  non- remote 
operations.  The  motor  is  attached  to  its  associated  firing  hardware  with 
cam-action  containter  latches.  The  motor  dimensions  are  0.750-in.  1,0, 

X  0,950-in.  O.  D.  with  the  same  lengths  as  for  the  threaded  motors. 
Shoulders  are  machined  on  the  outer  surface  near  the  motor  ends  to 
facilitate  casting-equipment  disassembly.  Motors  are  currently  fabricated 
from  both  stainless  steel,  and  phenolic-glass-fiber  tubes,2 

'  Detailed  shop  drawings  of  all  hardware  and  components  are  available 
upon  request. 

zTbe  glass-fiber  motors  were  purchased  from  Dynetics,  Inc.,  34  Crestview 
Road,  Mountain  Lakes,  New  Jersey. 
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The  firing  head  and  nozzle  each  contain  a  counterbore 
with  a  radial  groove  for  O-ring  sealing  on  the  external  motor  surface.  The 
nozzle  throat  inserts  are  graphite.  Anothe.r  O-ring  seals  the  firing  head 
to  the  pressure  transducer. 

Two  Camloc  light-duty  container  latches1  are  bolted  to 
the  nozzle  and  are  used  to  clamp  the  motor  assembly  together  during 
firing  (Fig.  2).  Extensions  on  the  Camloc  latches  permit  the  use  of  motors 
having  different  lengths,  and  reinforcing  ribs  welded  to  the  hooke  permit 
operation  at  a  maximum  motor-chamber  pressure  of  2000  psi. 

A  special  nozzle-closure  assembly  provides  reproducible 
blow-out  pressures.  The  design  includes  the  use  of  closure  discs  having 
thicknesses  from  0.001  to  0.050  in.,  a  metal  blow-out  slug  which  must 
shear  a  plug  from  the  cLosure  before  escaping,  and  a  closure  retainer 
which  screws  onto  the  nozzle  exit.  Blow-out  pressures  are  controlled  by 
the  selection  of  closure-disc  material  and  thickness. 

Stacking  sleeves  are  used  with  standard  motors  to  obtain 
additional  lengths  as  required,  but  a  special-length  clamp  extension  is 
required  for  each  combination.  A  dummy  load  cell  machined  from  stainless 
steel  is  used  with  a  commercial  pressure  transducer  when  firing  new 
propellants  which  have  unknown  sensitivity  characteristics. 

2.2.3  Clamp-Type  ZC1.5  -  faiotpy _ 

The  2C1.5  -  4-c.Iamp-type  motor  (Fig.  3). was  designed 
primarily  for  remote  operations,  but  its  ease  of  assembly  and  good  grain- 
dimension  reproducibility  makes  it  equally  useful  for  routine  operations. 
When  new  materials  having  unknown  sensitivity  characteristics  are  being 
tested  and  remote  operation  is  required,  a  motor2  made  of  filament- 
wound  epoxy— glass  tubing  is  used  to  minimize  potential  fragment  hazards 
from  malfunctions.  These  motors  are  wound  and  machined  to  2.000-in, 

I.  D.  X  2.258-in.  O.  D,  X  4-in.  long  with  a  shoulder  near  une  end  to  hold 

^Part  No.  51L-1-1AA,  Camloc  Fastener  Corporation,  Puramus,  New  Jersey. 
Purchased  from  Lamtex  Industries,  Inc.  ,  of  Farrningdale,  Long  Island, 

New  York, 
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Fig  3  Clamp-type  PCI  5- <4  motor  firing  assembly 
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the  motor  during  casting-equipment  disassembly.  The  outer  surface  of 
the  motor  is  used  as  an  O-ring  seal  and  casting-fixture-alignment  surface. 

For  non-remote  operation  a  motor  made  of  2.000-in,- I.  D,  stainless  steel 
tubing  having  dimensions  similar  to  those  of  the  glass  fiber  motor  may 
be  used. 

The  firing  head  contains  a  counterbore  with  an  O-ring 
groove  for  inserting  one  end  of  the  motor,  two  threaded  ports  for  installing 
commercial  pressure  transducers,  a  threaded  boss  for  mounting  to  a 
commercial  load  cell,  and  two  Camloc  heavy-duty  container  latch1  assemblies 
for  clamping  the  motor  firing  asserribly  together.  Two  mating  Camloc 
hooks2  are  welded  to  the  nozzle  housing  which  also  contains  a  counterbore 
with  an  O-ring  groove  for  inserting  the  other  end  of  the  motor.  Both 
firing  head  and  nozzle  housing  are  made  of  mild  steel.  Steel  nozzle 
inserts  containing  press-fitted  graphite  throats  are  sealed  in  the  nozzle 
housing  with  O- rings  and  are  retained  by  a  threaded  steel  ring. 

The  Camloc  latch  as  received  from  the  manufacturer 
will  function  at  a  maximum  chamber  pressure  of  1700  psi,  but  the. 
maximum  chamber  pressure  can  be  increased  to  2000  psi  by  replacing 
the  latch  shear  pins  with  hardened  steel  bolts.  Latch  extension  arms 
serve  as  turnbucklea  for  adjusting  the  latch  length  and  are  fabricated  in 
length  increments  required  for  use  with  4-in.  ,  8-in,  ,  12-in.  ,  and  16-in. 
motor  combinations.  All  motors  are  cast  in  4-in.  lengths;  the  additional 
lengths,  which  are  used  only  for  non-remote  firings,  are  obtained  by 
coupling  the  motors  end  to  end  through  the  use  of  O-ring-sealed  short 
sleeves  which  provide  a  '/16-in.  gap  between  motor  sections.  The  use  of 
turnbuckle-type  extensions  requires  fabrication  of  latch  lugs  with  left-  hand 
threads.  Both  the  lugs  and  extensions  are  fabricated  from  Type  4130 
steel  and  heat  treated. 


'Camloc  Part  No.  37L11-1B. 
2Cainloc  Part  No.  37L15-1A. 


A  special  firing  head  and  nozzle  housing,  with  the  latch 
assembly  replaced  by  y2-in.  steel  tie  rods,  are  used  for  non-remote 
assembly  of  motors  when  chamber  pressures  from  2000  to  3000  psi  are 
anticipated.  A  dummy  Load  cell  is  used  when  firing  novel  propellants 
which  have  unknown  characteristics. 

2.3  Remote  Casting  Equipment  and  Facilities 

2.3.1  .75C.50  Motor  Casting  Fixtures 

Casting  fixtures  (Fig.  4)  for  both  the  threaded-  and  clamp 

type  .75C.50  rhotors  are  assembled  and  disassembled  using  a  screw- 

clamp  fixture.  Two  stationary  yokes  at  the  top  of  the  clamp  hold  the 

motor  shoulder  and  casting  head,  and  a  screw-actuated  yoke  holds  the 

one-piece  casting  base  and  mandrel.  Mandrel  alignment  is  provided  by 

a  Teflon  gasket  at  the  mandrel  base,  and  by  a  spacer  at  the  top,  which 
-  •  - i  •• 

is  inserted  temporarily  during  assembly  and  removed  immediately  prior 

to  casting.  The  clamp,  with  the  screw  manually  driven  by  a  hand  crank 

or  remotely  driven  by  a  geared  air  motor,  disassembles  the  motor  and 

casting  fixtures.  A  Teflon  plug  fits  into  the  upper  end  of  a  0,50-in.-O.  D. 

hollow  mandrel  during  the  filling  operation.  When  filling  is  completed, 

the  plug  is  removed  and  a  Teflon  overfill  plug  is  inserted  which  forces 

the  excess  propellant  into  the  hollow  portion  of  the  mandrel.  The  overfill 

plug  in  conjunction  with  a  gasket  at  the  base  of  the  mandrel  provides  a 

grain  recess  at  each  end  of  the  motor. 

2.3.2  2C1.5-4  Motor  Ca sting  Fixtures 

Casting  fixtures  for  the  2C1.5  motor  (Fig.  5)  consist  of 
a  casting  head,  a  one-piece  casting  base  and  mandrel,  a  break-off 
ring,  and  two  tie  rods.  Small  wooden  cross  pins  are  used  on  the  bottom 
ends  of  the  tie  rods  to  act  as  shear  pins  during  remote  disassembly.  A 
shoulder  on  the  base  of  the  mandrel,  ar.d  the  break-away  ring  provide  a 
l/l(,-in.  grain  recess  at  each  end  of  the  motor.  A  remote  casting-fixture 
disassembly  device  (Fig.  6)  actuated  by  a  double-acting  air-draulic 
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Fig.  5  Casting  equipment  for  2 C  1,5-4  motors. 

cylinder  is  used  to  disassemble  the  motor.  The  casting  fixture  is 
set  on  the  jig  in  an  inverted  position.  Aa  the  moveable  portion  comes 
down,  two  lingers  press  against  the  tie  rods,  removing  them,  the  casting 
head,  and  the  break-off  ring.  The  puller  jaw  engages  a  button  on  the  end 
of  the  mandrel,  and  the  mandrel  is  extracted  on  the  upward  stroke. 

2.3.3  Remote  Casting  Facilities 

luree  facilities,  differing  primarily  in  size  and  propel¬ 
lant-handling  capacities,  were  designed  and  constructed  for  mixing  and 
casting  hazardous  propellants.  Two  small  facilities  were  built  to  handle 
up  to  2  50  g  of  propellant.  Operations  are  carried  out  behind  a  semicircular 
steel  shield  (Figs.  7  and  8)  having  a  radius  of  three  feet  with  two  AMF 
Mini-Manip  manipulators1;  a  rectangular  Plexiglas2  sight  port  provides 
for  direct  viewing.  The  air  within  each  shield  may  be  maintained  at  a 
dew  point  of  -50°C,  thus  allowing  the  processing  of  hygroscopic  materials. 

'Part  No,  89-63-R1000,  AMF  Atomics,  A  Division  of  American  Machine  & 
Foundry  Company,  140  Greenwich  Ave.  ,  Greenwich,  Conn, 
zRohm  Sj  Haas  Company,  Philadelphia ,  Pennsylvania. 
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One  shield  was  designed  for  casting  .7  5C.50  motors  and  the  other  was 
designed  for  casting  both  .75C.50  and  2C1.5  motors. 

The  third  facility  (Figs.  9  and  10)  was  designed  to 
handle  up  to  2  lbs  of  propellant.  Mixing,  casting,  curing,  and  motor 
disassembly  are  carried  out  remotely  in  a  concrete  cell.  Transfer 
operations  are  handled  with  two  AhiF  hiodel  8  extended  — reach  standard- 
duty  manipulators  while  processing  equipment  is  operated  with  compressed 
air  or  hydraulically.  A  sight-port  containing  a  3 -in.  block  of  Plexiglas 
is  used  for  direct  viewing,  and  close-up  views  of  equipment  and  controls 
are  provided  with  a  closed-circuit  television  camera. 
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2 . 4  Remote  Motor  Assembly  and  Firing  Facility 

The  remote  motor  assembly  and  firing  facility  was 
designed  for  installation  within  an  existing  static-firing-range  bunker. 

The  facility  provides  for  completely  remote  assembly,  firing,  or 
disassembly  of  .75C.50  and  2C1.5-4  clamp-type  motors  without  personnel 
exposure. 

A  thrust  plate  for  firing  motors  in  a  vertical  position  is 
mounted  at  table-top  level  on  a  reinforced  concrete  block.  A  1-in.  -thick 
steel  shield  is  set  in  front  of  the  stand  and  a  heavy  wire-mesh  screen 
extends  from  the  top  of  the  shield  to  the  bunker  ceiling  to  protect  the 
operator  while  the  motor  is  being  set  on  the  firing  stand.  Transfer 
operations  within  the  test  cell  are  handled  with  a  Mini-Manip  manipulator 
which  is  set  in  the  upper  right-hand  corner  of  the  shield  while  motoirs. 
are  prepared  for  firing.  A'  14-in, -wide  X  22-in. -high  opening  is  provided 
in  the  shield  for  a  movable  Plexiglas  sight-port  which  also  serves  as  an 
access  door.  The  door  is  air-powered  and  is  remotely  actuated. 

Motors  are  brought  to  the  stand  for  firing  in  an  armoured 
transporter.  An  air-cylinder-actuated  retracting  cup  on  the  stand  is 
automatically  actuated  by  door  movements.  When  the  door  is  opened, 
the  cup  swings  toward  the  door  opening.  The  motor  is  transferred  from 
the  transporter  to  the  cup  with  a  Mini-Manip  on  the  transporter.  The 
door  is  closed,  and  the  cup  and  motor  are  retracted  inside  the  shield. 

Fig.  11  shows  an  outside  view  of  the  shield  with  the  door  open  and  the 
manipulator  installed.  With  the  shield  door  closed,  the  operator 
assembles  the  motor  remotely,  removcs-the  manipulator,  and  leaves 
the  bay  before  the  motor  is  fired. 

The  ,7  5C.  50-motor  assembly  and  disassembly  fixtures 
use  a  pivoted  motor-clamping  yoke  and  a  pivoted  unlatching  hook  in 
conjunction  with  the  manipulator.  Fixtures  for  the  remote  assembly  and 
disassembly  of  20.5  -  4  motors  are  clamped  and  unclamped  by  ropes 
extending  through  the  front  of  the.  shield.  Figs.  12  and  13  illustrate  some 
steps  in  the  assembly  of  .75C.30  and  20. 5  damp-type  motors. 


Best  Available  Copy 


(a)  GAUGES,  FIRING  HEAD  (B)  MOTOR  INSERTED  INTO 

AND  IGNITER  IN  PLACE  FIRING  HEAD 


(c)  NOZZLE  ASSEMBLY  (d)  COMPLETE  ASSEMBLY, 

READY  TO  FIRE 


Fig.  12.  Steps  in  Ihe  assembly  of  the  75C.5C  micro-motor  for  firing  in  the 
static  lest  stand. 


Best  Available  Copy 


(a)  GAUGES  AND  FIRING 
HEAD  IN  PLACE 


(B) MOTOR  CASE  IN  PLACE 
NOZZLE  BEING  TRANS¬ 
FERRED. 


(d)  COMPLETE  ASSEMBLY 
READY  TO  FIRE. 


the  20. 5-4  motor  for  firincj  on  the 
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The  facility  was  tested  by  detonating  250  g  of  C-4 
explosive  confined  in  a  2C1, 5-4 fiberglass  motor  located  at  the  position  of 
normal  motor  firing.  The  sight-port  withstood  a  separate  test  with 
250  g  of  C-4  in  a  fiber glas s  motor  located  3  in.  from  the  inside  pane. 

2.5  Self-Propelled  Hazardous  Materials  Carrier 

A  powered  vehicle  (Fig.  14)  capable  of  safely  transporting 
1  lb  of  hazardous  material  from  the  remote-casting  facilities  to  the 
remote-firing  facility  without  exposure  of  personnel  was  designed, 
fabricated,  and  tested.  The  design  was  patterned  after  a  hazardous 
material  vehicle  designed  and  built  by  the  E.  I.  duPont  de  Nemours  Sc 
Company,  Inc.  1  The  vehicle  consists  of  a  warehouse-type  electric  truck 
chassis2  upon  which  is  mounted  a  30-in.  sphere  made  of  1.5-in.  steel. 

The  vehicle  is  operated  from  a  shielded  cab. 

The  sphere  is  equipped  with  hydraulically  actuated 
devices  which  open  an  inwardly  swinging  lid  and  extend  a  cup  into  which 
hazardous  material  is  placed.  Transfer  operations  are  handled  remotely 
with  a  modified  Mini-Manip  manipulator,  installed  in  the  front  of  the  cab. 
The  walls  of  the  cab  are  1  in,  thick  in  front  (the  surface  most  directly 
exposed  to  the  sphere)  and  l/z  in.  thick  on  the  sides  and  top.  Aluminum 
doors  in  the  rear  protect  the  operator  from  blast  overpressure  in  the 
event  of  an  explosion  outside  of  the  sphere. 

2 . 6  .75C.50  Motor  Shielded  Transport  Containers 

Two  types  of  portable  containers  are  used  to  transport 
.75C.50  motors:  a  thick-wall  closed  vessel,  and  an  open-end  pipe-section 
carrier  (Figs.  15  and  16). 


1 


'Explosives  Department,  Repauno  Development  Laboratory,  Gibbstown, 
New  Jersey. 

^Model  No.  RW-41136,  Automatic  Transportation  Company,  Chicago, 
Illinoi  s . 
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3,  Experimental  Techniques 

3 . 1  Motor  Preparation  and  Physical  Measurements 

3.1.1  Steel  Motor  Cleaning 

Steel  motor  cases  are  hydroblasted  thoroughly,  rinsed 
with  water  to  remove  all  grit,  dipped  in  acetone  to  remove  the  water, 
and  quickly  transferred  to  a  trichloroethylene-vapor  degreaser.  After 
remaining  in  the  degreaser  for  several  minutes,  the  motors  are  removed, 
capped  on  each  end  with  plastic  cap-plugs,  and  stored  at  a  temperature 
of +140*  F  until  required  for  use. 

3.1.2  Fiberglass  Motor  Cleaning 

Fiberglass  motors  are  hydroblasted  very  lightly  to 
remove  the  surface  glaze.  The  remaining  operations  are  the  same  as 
for  the  steel  motors,  except  that  time  in  the  degreaser  is  reduced  to 
about  30  seconds. 


Fiq  15  Deton.obie  sample  container 
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Fig. 16  Open-end  micro-motor  carrier . 
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3.1.3  Lining 

The  motor  liner  depends  on  the  type  of  propellant  to 
be  cast.  For  composite  double-base  plastisol  propellants,  the  cases 
are  sprayed  with  an  extremely  thin  film  of  Rohm  &  Haas  PL-1  liner 
(Table  I).  Many  experimental  propellants  will  bond  adequately 
to  the  clean  steel  or  fiberglass  motor  cases,  and  no  liner  is  required. 

Any  liner  must  be  thin  and  uniform  to  avoid  errors  in  burning  rate 
measurements,  as  well  as  efficiency  losses  associated  with  long  tail-offs. 

Table  I 

Composition  of  PL-1  Liner 
Constituent  %  by  Weight 


Cellulose  Acetate  (Powdered) 

4.57 

Tri phenyl  Phosphate  Flakes 

2.88 

Santicizer  M-17*5 

1.74 

Red  Lead  (Pb304) 

.02 

Tolylene  Diisocyanate  (Hylene  TM) 

.23 

Acetone  Tolylene 

46. 3y 

Methyl  Cellos  olve  Acetate^ 

44.17 

j^Type  384,  Tennessee  Eastman  Co.  ,  Kingsport,  Tenn. 

Monsanto  Chemical  Co.,  St.  Louis,  Mo, 
c 

^E.  I.  duPont  de  Nemours  Co.  ,  Wilmington,  Del. 

'  Ethylene  Glycol  Monomethyl  Ether  Acetate,  Union  Carbide 

Chemical  Co.,  New  York,  N.  Y. 

3.1,4  Casting 

While  propellant  charges  for  rnicro-motors  could  be 
formed  in  a  number  of  ways,  all  motors  tested  by  this  Division  have  been 
loaded  by  casting  the  propellant  directly  into  the  motor  case  and  curing 
it  in  place. 

Casting  techniques  vary  widely  depending  on  the  specific 
propellant  composition.  For  example,  compositions  which  are  very  fluid 
may  be  cast  easily,  while  extremely  viscous  compositions  present 


special  problems.  Specific  casting  techniques  are  described  in  two  aeries 
of  reports  of  this  Division:  Quart erly  Progress  Reports  on  Chemical 
and  Propellant  Processing,  and  Quarterly  Progress  Reports  on  ARPA 
Projects. 

3.1.5  Motor  Cleanup 

A  certain  amount  of  motor  cleanup  is  usually  required 
after  the  grain  has  been  cast.  All  external  motor  surfaces,  especially 
those  which  mate  with  the  nozzle  and  head  closure,  must  be  thoroughly 
cleaned.  A  solvent  such  as  acetone  is  usually  used  for  this  purpose. 

3.1.6  Grain  Trimming  and  Inspection 

The  goal  during  casting  is  the  production  of  a  perfect, 
cast-to-length  grain.  However,  it  is  sometimes  necessary  to  trim  the 
grain  to  remove  end  flaws.  In  this  event  the  motor  is  clamped  in  a  vise 
and  a  small  90*  knife  is  used  to  trim  the  grain  ends,  using  a  special 
trimming  jig  as  the  knife  guide.  After  trimming  and  cleaning,  the  motor 
is  visually  inspected  for  propellant  voids,  crackB,  case-bond  failures  or 
other  irregularities. 

3.1.7  Weighing 

Motor  cases  are  weighed  after  lining,  before  firing 
(after  trimming),  and  after  firing.  The  .75C.50  motors  are  weighed  on  a 
200-gram  direct-reading  Mettle r  electric  analytical  balance  which  can 
be  read  to  0.1  milligram.  The  ZC1.5-4  motors  are  weighed  on  a  similar 
800-gram  Mettle r  balance  which  can  be  read  to  0.01  gram.  Due  to  the 
small  size  of  the  motors,  weights  must  be  determined  with  extreme 
accuracy  to  provide  accurate  specific-impulse  measurements.  Conse¬ 
quently,  particular  care  must  be  taken  to  insure  that  the  motor  case 
before  and  after  firing  is  unaltered  with  the  exception  of  propellant  and 
liner  consumed. 
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3.1.8  Grain  and  Nozzle  Measurements 

Grain  length  is  measured  to  0.01  inch  with  a  standard 
steel  rule.  Grain  web  thickness  is  measured  to  0.001  inch  from  measure¬ 
ments  of  the  grain  I.  D.  before  firing  and  case  I.  D.  after  firing;  the 
motor  is  not  cleaned  before  the  latter  measurement.  Nozzle-throat 
diameters  are  measured  to  0.001  inch  before  firing,  after  firing  on  top 
of  any  material  deposited  in  the  nozzle,  and  after  removal  of  deposits 
in  the  nozzle.  Small  expanding -ball-type  hole  gauges  and  standard 
micrometers  are  used  for  throat-diameter  measurements.  Nozzle 
exit  diameter  is  measured  to  0.001  inch  using  an  inside -reading  caliper 
micrometer. 

Extreme  care  must  be  taken  to  determine  such  physically 
measured  quantities  as  web  thickness,  nozzle  throat  and  exit  diameter, 
and  charge  weight,  in  order  to  obtain,  acceptable  ballistic  data.  Due  to 
the.  small  magnitude  of  these  quantities,  allowable  tolerances  on  measure¬ 
ments  are  inherently  small  and  difficult  to  maintain. 

3.  2  Motor  Assembly 

3.2.1  Threaded  .75C.50  Motors 

For  threaded  .75C.50  motors,  the  special  pressure 
gauge  serves  as  the  motor  head-closure  and  contains  the  squib  wire 
pass-through  port.  The  igniter  is  installed  by  passing  the  squib  wires 
down  into  the  center  of  the  gauge  and  out  through  the  igniter  port,  threading 
the  wires  through  the  close -tolerance  holes  in  the  Teflon  compression 
gland,  pulling  the  igniter  well  down  to  position  it  at  the  motor  head  end, 
and  inserting  and  tightening  the  ignition  port  screw  to  provide  the 
pressure  seal.  A  disc-type  nozzle  closure  is  installed  at  the  nozzle 
exit  using  a  closure-reitaining  ring  which  threads  onto  the  nozzle  aft  end. 
The  pressure  gauge  and  nozzle  are  threaded  onto  the  ends  of  the  motor 
to  complete  the  loading  operations.  A  sectioned  view  of  the  loaded  motor  ' 
is  shown  in  Fig.  1 . 
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3.2.2  Clamp-Type  .75C.50  Motors 

For  the  clamp-type  .75C.50  motors,  the  igniter  is 
installed  in  the  motor  firing  head  in  a  manner  completely  analogous 
to  that  used  for  igniter  installation  into  the  threaded  motor  pressure 
gauge.  The  disc -type  nozzle  closure  is  similarly  retained  at  the 
nozzle  exit  by  a  threaded  retainer  ring.  Two  types  of  nozzle -closure- 
retainer  ring  assemblies  are  available:  one  clamps  the  closure  to  the 
exit  (as  for  the  threaded  motor),  and  one  utilizes  a  coin-type  slug  to 
rupture  the  closure,  allowing  better  control  of  closure  rupture  pressure. 
The  reinforced  draw  hooks  are  attached  to  the  nozzle  clamp  arms  (with 
or  without  extensions,  depending  on  motor  length)  and  adjusted  in  length 
for  proper  assembly  tension.  Final  assembly  is  made  by  inserting  the 
firing  head  into  the  special  clamp-type  pressure  gauge,  inserting  one 
end  of  the  motor  into  the  firing  head,  pushing  the  nozzle  down  over  the 
motor  end,  placing  the  draw-hook  feet  under  the  rim  of  the  pressure  gauge 
(clamp  levers  in  a  raised  position),  and  seating  the  clamp  lovers.  The 
0-rings  and  the  outer  motor  ends -should  be  lightly  coated  with  a  light 
vacuum  grease  and  care  must  be  taken  not  to  pinch  the  nozzle  and 
firing-head  O-rings  during  this  operation.  Since  the  pressure  seal  is 
made  by  the  O-rings  on  the  motor  case  O.  D. ,  clamp-arm  tension  is 
not  critical  and  need  only  be  sufficient  to  provide  a  snug  assembly,  A 
sectioned  view  of  the  assembled  motor  is  shown  in  Fig.  2. 

3.2.3  Clamp-Type  2C1.5-4  Motors 

For  the  clamp-type  2C1.5-4  motors,  the  igniter  and 
closure  are  installed  in  the  nozzle  insert,  which  is  then  pushed  into  the 
nozzle  shell  and  secured  by  threading  down  the  nozzle  retaining  hut. 

The  disc -type,  bevelled  closure  is  glued  in  the  nozzle  convergent  section. 
The  squib-wire  length  is  adjusted  so  that  the  igniter  will  be  positioned  at 
the  motor  head-end  when  the  motor  is  assembled.  The  outside  ends  of 
the  motor  and  the  nozzle  and  firing-head  O-rings  are  lightly  coated  with 
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light  vacuum  grease.  Final  assembly  is  made  by  inserting  one  end  of 
the  motor  into  the  firing  head,  pushing  the  nozzle  housing  down  over  the 
motor  end,  placing  the  clamp  lugs  in  the  nozzle  housing  ears,  and 
seating  the  clamps.  As  for  the  clamp-type  .75C.50  motors,  clamp  arm 
tension  is  not  critical  and  need  only  to  be  sufficient  to  provide  a  9nug 
assembly.  A  section  view  of  the  assembled  motor  is  shown  in  Fig.  3. 

3. 3  Ignition 

3.  3.  1  Ignition  of  „75C.5Q  Motors 
3.3.  1.1  General 

One  of  the  major  problem  areas  in  the  .75C.50  motor 
program  has  been  the  development  of  a  satisfactory  ignition  system. 

Many  different  types  of  ignition  systems  have  been  tested,  utilizing  a 
variety  of  squibs,  igniter  materials,  igniter  weights,  igniter  positions, 
nozzle  closures,  and  closure  positions.  Conclusions  drawn  as  a  result 
of  these  tests  are  summarized  below. 


3.  3. 1.2  Nozzle  Closures 

A  nozzle  closure  is  essential  to  reduce  ignition  delay 
time  and  rise  time,  and  to  eliminate  misfires  and  hangfires.  For  best 
results,  the  closure  should  rupture  at  pressures  between  80%  and  120% 
of  the  motor  equilibrium  pressure.  Nozzle  closures  placed  in  the  nozzle 
convergent  section  have  usually  caused  severe  ignition-pressure  spikes; 
much  better  results  are  obtained  with  a  disc -type  nozzle  closure  clamped 
to  the  nozzle  exit. 

Disc -type  nozzle  closures  should  fail  by  shearing 
cleanly  at  a  diameter  larger  than  the  nozzle  exit  diameter.  If  the  closure 
docs  not  shear  cleanly  rupture  pressures  will  be  variable  and  thrust 
measurements  will  be  affected  by  gas-flow  disturbances  produced  by 
protrusion  of  jagged  edges  of  the  closure  into  the  exit  gas  stream.  To 
■insure  clean  shearing  of  the  nozzle-  closure,  a  special  closure-retainer 
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is  available  which  utilizes  a  coin-type  slug  to  forcibly  rupture  the  closure 
{Fig.  2).  This  assembly  also  provides  much  better  precision  of  closure 
blowout  pressure  when  metal  closures  are  used. 

Nominal  nozzle-closure  blowout  pressures  for  regular 
variable -diameter  retainer  rings  are  given  in  Table  II  and  nominal 
nozzle -closure  blowout  pressures  for  the  special  plug -type  closure- 
retainer  assembly  are  given  in  Table  III. 


Table  II 

Nozzle -Cloauro  Blowout  Pressures  for  Regular  Retainers 


Closure  Material 
and  Thickness 

Closure  Blowout 

pr„  ...... „a  „ 

0.625"  Din 
Retainer 

rector  0.75 n  Dia motor 
Retainer 

0,875M  Diameter 
Retainer 

1.00"  Diameter 
Retainer 

CAb-.0H)*< 

480 

400 

340 

300 

0  A  -.015" 

700 

590 

500 

440 

CA  -.020" 

920 

780 

670 

580 

CA  -.025" 

1150 

970 

830 

720 

CA  -.030" 

1380 

1160 

1000 

870 

CA  -.035" 

1620 

1350 

1 160 

1010 

CA  -.040" 

1850 

1540 

1320 

1160 

CA  -.045* 

2080 

1730 

1490 

1300 

CA  -.050" 

2310 

1920 

1650 

1450 

CA  -.060" 

2770 

2310 

1990 

1730 

^Itr  3td.  deviation  c  200  pula. 
Cellulose  acetate  sheet. 
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Nozzle-Closure  Blowout  Pressures  for  Plug-Type  Assembly 


Material 

Thickness,  in 

Blowout  Pressure,  paia 

1(T  ,  psia 

Aluminum3, 

0.008 

480 

7 

Aluminum 

0.010 

600 

27 

Aluminum 

0.012 

700 

24 

Aluminum 

0.016 

930 

33 

Aluminum 

0.020 

1170 

34 

Aluminum 

0.025 

« 

1400 

26 

CAb 

0.020 

750 

202 

CA 

0.025 

1060 

139 

CA 

0.030 

1360 

115 

CA 

0.035 

1540 

247 

c 

Brass 

0.005 

1060 

60 

100-0  aluminum  sheet. 
Cellulose  acetate  sheet. 

c 

Sheet  shim  stock. 


.  .  3.  3.  1. 3  Igniters 

The  optimum  igniter  location  is  at  the  head  end  of  the 
motor  with  the  squib  wires  passing  through  the  firing  head.  If  the  Igniter 
is  toward  the  aft  end,  the  squib  head  may  be  torn  from  the  squib  wires 
and  plug  the  nozzle.  Satisfactory  ignition  can  usually  be  accomplished 
with  an  igniter  positioned  in  the  nozzle  expansion  cone,  but  the  necessary 
piercing  of  the  nozzle  closure  to  permit  squib-wire  passage  results  in  ■ 
ragged  closure  rupturing  and  erratic  rupture  pressures. 


. . .  wMwn 
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For  propellants  which  are  relatively  easy  to  ignite,  the 
optimum  igniter  consists  of  an  M-l  squib  for  the  1.5-inch  motor,  or  a 
squib  with  additional  RIP-11  or  RHim-1  (Table  IV)  for  the  2,5  and  3.5 
inch  long  motors.  The  larger  igniters  are  made  by  pouring  the  required 
amount  of  igniter  powder  into  the  small  half  of  a  number  0  gelatin  caps\ile, 
inserting  the  squib  into  the  capsule,  and  cementing  the  squib  and  capsule  to¬ 
gether  with, DiicO3 cement.  The  resultant  igniters  are  of  the  fast-acting,  flash 

Tabic  IV 

Composition  of  RHjm-1  Igniter  Powder 


Constituent _ Weight  per  cent 


Metal  Magnesium  (55  —  1  Opjo.) 

6o 

KCIO*  (~105p) 

25 

Ba(NQ3)2  (<  149(a) 

15 

For  propellants  which  cannot  be  satisfactorily  ignited 
with  the  M-l  squib -RHim  -1  system,  a  dipped  Atlas  Match  igniter  is 
generally  satisfactory.  This  igniter  consists  of  a  plain  Atlas  Match 
squib  which  is  repeatedly  dipped  into  a  slurry  of  RHim-1  and  a  solution 
of  5%  poly(isobutylene)  in  n-hexane  until  a  sufficient  quantity  of  the  mixture 
adheres  to  the  squib  head.  The  solution  is  then  allowed  to  evaporate, 
leaving  a  coating  of  RHim-1  on  the  squib  head.  The  resultant  igniters 
are  of  the  soft,  slow -burning  type,  and  eject  burning  particles. 


3,3.  1,4  Effect  of  Ignition  on  Specific- Impulse  Measurements 


Corrected  specific  impulse  using  the  .75C.50  motors  is 
markedly  dependent  upon  ignition  delay  time,  ignition  weight,  and  ignition 
rise  time.  Long  ignition-delaytimes  invariably  result  in  low  values  of 


^IF-l  is  a  Picatinny  Arsenal  product  consisting  of  68,2%  barium  nitrate 
and  31.8%  magnesium. 

"E,  I.  DuP'int  Hp  Nornoii'-s  &  Co.  ,  Wilmington.  Delaware. 
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specific  impulse  due  to  propellant  cook-off  and  increased  heat  losses. 

On  the  other  hand,  the  measured  specific  impulse  will  be  too  high  if  a 
relatively  large  weight  of  igniter  material  is  U9ed,  and  no  correction  is 
made  for  the  contribution  of  the  igniter.  Measurements  of  specific  impulse 
of  a  typical  plastisol  propellant  as  a  function  of  ignition  delay  time  and 
igniter  weights  demonstrated  the  need  for  short  ignition-delay  times  with 
minimum  quantities  of  igniter  material  (Fig,  17). 

The  policy  of  this  Division  has  been  to  use  igniters  which 
weigh  no  more  than  one -half  per  cent  of  the  charge  weight,  and  to  make  no 
igniter  correction  to  specific  impulse  although  the  igniter  adds  about  %% 
to  total  impulse.  This  increase  is  very  closely  matched  by  an  apparent 
decrease  in  impulse  due  to  the  fact  that  motors  are  fired  vertically  and 
the  loss  of  weight  results  in  a  base-line  shift  of  the  thrust  cell.  The  M-l 
squib  used  to  ignite  the  1.5 -inch  motor  contains  about  50  milligrams  of 
combustible  material,  which  is  one-half  per  cent  of  the  nominal  charge 
weight,  .'or  2.5-  and  3,5-inch  motors,  RHim-1  igniter  powder  is  added  to 
make  the  total  igniter  weight,  including  the  M-l  squib  weight,  equal  one- 
half  per  cent  of  the  charge  weight. 


200  400  600  600  1000  1200  1400  I6C0 

IGNITION  DELAY  TIME  -  msec 

Fig  17  Ft  feci  of  ignition  delay  and  igniter  •vsiqKi  on  the  specific  impulse 
meosuied  m  the  75 C  50  micro  motor 


-33- 


If  the  weight  of  the  igniter  necessary  to  achieve  ignition 
is  greater  than  one-half  per  cent  of  the  charge  weight,  an  igniter  correction 
should  be  applied  to  the  measured  specific  impulse.  The  method  used  is 
discussed  in  Section  3.6. 

3.3.  1.5  Effect  of  Ignition  on  P-K-r  Measurements 

Ignition  requirements  for  P-K-r  measurements  differ 
somewhat  from  those  for  specific-impulse  measurements.  The  primary 
objective  is  a  pressure  trace  in  which  this  pressure  builds  up  rapidly,  end- 
of -burning  is  sharp,  and  pressure  decays  rapidly.  To  obtain  this,  an 
M-l  squib  with  additional  RHim -powder  is  used.  Since  the  igniter-correction 
problem  which  must  be  considered  in  specific-impulse  measurements  is 
not  critical  for  P-K-r  measurement*,  any  amount  of  additional  RHim-1  required 
may  be  used  to  achieve  the  desired  trace  shape.  In  most  cases  .05  to 
.15  grama  of  RHim-1  is  sufficient  for  ignition  of  the  1.5 -inch  motor.  Clean 
shearing  of  nozzle  closures  is  not  critical. 

3.3.2  Ignition  of  2 C  1.5 -4  Motors 

3. 3. 2.1  General 

Two-inch-diameter  motors  are  easier  to  ignite  than  the 
.75 -inch  because  as  much  as  0.7  gram  (%%  of  the  charge  weight)  of 
igniter  material  can  be  used  for  specific-impulse  measurements.  Larger 
igniters  may  be  used  for  P-K-r  measurements. 

3.3.  2.2  Nozzle  Closures 

Closures  which  rupture  at  pressures  from  20  to  100  per  cent 
of  motor  operating  pressure  are  generally  acceptable  for  2 -inch  motors. 

Plastic  nozzle  closures  are  usually  cemented  into  the  nozzle  convergent 
section  with  an  adhesive  such  as  Pliobond^1.  An  injection-molding  process 
is  used  to  fabricate  these  closures  (Fig.  18)  in  three  diameters  and  a  wide 
range  of  thicknesses.  Closures  having  a  thickness  between  0.035  and 
0.060  inch  generally  provide  excellent  ignition  for  firings  at  a  nominal 
pressure  of  1000  psia. 

'The  Goodyear  Tiro  &  Rubber  Co.,  Akron  1  6,  Ohio. 
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10 -I  10-2  10-3 

Fig, IQ  Nozzle  closures  for  2C.I.5-4  moiors. 


3.  3.  2.3  igniters 

Two  types  of  igniters  have  been  used  for  2-inch  motors: 
a  fast-acting  bag  type  and  a  Blow-burning  jelly-roll  type.  Both  utilize 
the  Atlas  Match  squib. 

The  bag-type  igniter  is  made  by  enclosing  the  igniter 
powder  (either  RIP  or  RHim)  and  the  Atlas  Match  squib  head  in  a  2 -mil 
polyethylene  bag.  The  bag,  cut  from  sheet  stock,  is  fastened  around  the 
squib  wires  with  a  rubber  band. 

The  jelly-roll  igniter  is  made  by  rolling  up  a  narrow 
rectangular  strip  of  the  previously  prepared  jelly-roll  material  around 
the  squib  head;  this  is  then  wrapped  tightly  with  one  layer  of  2 -mil  poly¬ 
ethylene  and  secured  with  rubber  bands.  The  finished  product  resembles 
a  jelly-roll,  with  the  squib  head  in  the  center.  The  jelly-roll  material 
iB  prepared  by  coating  a  sheet  of  cheesecloth  to  a  depth  of  about  Vi^-inch 
with  a  slurry  of  RHim-1  and  a  solution  of  5%  poly(isobutylene)  in  n-hexane, 

3, 4  Insulation  of  .75C.5Q  Motors 

Quite  accurate  P-K-r  and  temperature-coefficient  data 
may  be  obtained  from  .75C.50  motors  if  the  propellant  temperature  is 
accurately  known.  However,  because  of  their  small  mass,  these  motors 
can  undergo  rapid  temperature  changes.  Consequently,  if  firing  temperatures 
other  than  ambient  arc  required,  the  motors  must  be  well-insulated  to  prevent 
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prohibitive  propellant-temperature  changes  from  occurring  during  the 
time  interval  between  removal  from  the  conditioning  box  and  firing. 
Temperature  measurements  were  made  by  inserting  a  thermocouple  in 
the  propellant  of  a  .75C. 50-1.5  motor.  The  motor  was  cooled  to  -40*F 
and  removed  from  the  cold  box.  After  four  minutes,  the  average  time 
from  box  to  firing,  the  temperature  had  risen  40*F  with  an  initial 
temperature  differential  of  125*F.(Fig.  19). 

Very  good  temperature  control  was  achieved  by  insulating 
the  loaded-motor  assembly,  including  the  pressure  gauge,  in  a  close- 
fitting  rectangular  box  made  of  1 -inch-thick  polystyrene  foam  (Fig.  20), 
The  box  is  open  on  the  pressure  gauge  end,  and  is  closed  on  the  nozzle 
end.  The  insulating  box  is  left  in  place  during  firing,  and  upon  ignition 
a  press-fit  plug  is  blown  out  of  the  closed  end  of  the  box. 


Fig  1 0  Chanqe  in  propellant  temperature  in  75C.5H.5  micro-motors  as  a 
function  of  time  after  removal  from  o  cold  box  at  -40  °F. 


**  m 1 1  u'Jim  ajiiii—juj»m,uiy«j  uijuum— i 
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Rg.20  Insulating  box  for  .75C.50-I.5  motor.  Assembled  box  and  motor 
at  left.  Cutaway  view  at  center  and  right. 


The  complete  insulated  assembly  is  conditioned  at  the 
required  temperature.  Several  of  these  assemblies  are  then  placed  in 
a  larger  insulated  box,  which  has  also  been  conditioned  to  the  motor 
temperature,  for  transportation  to  the  firing  range.  The  insulated  motor 
assemblies  are  removed  one  at  a  time  for  firing;  a  firing  can  usually 
be  accomplished  within  two  minutes  of  removal  from  the  insulated  trans¬ 
port  box. 

Temperature -time  measurements  made  with  a  thermocouple 
embedded  in  the  propellant  showed  that  the  insulation  reduced  heat  transfer 
sufficiently  that  accurate  estimates  of  the  propellant  temperature  at  the 
time  of  firing  could  be  made  (Fig.  19). 

Micro-motors  have  been  successfully  fired  at  extreme 
temperatures  of  -40*F  and  +140*F.  Since  the  pressure  gauges  are  con¬ 
ditioned  vdth  the  motors  they  must  be  calibrated  at  these  temperatures. 


Best  Available  Copy 


Since  the  standard  Celvacene1  light  vacuum  grease  used  to  pack  the  gauges 
is  satisfactory  at  +140*F  but  freezes  at  -40*F,  Dow-Corndng  High-Vacuum 
Grease  diluted  with  Dow  Silicone  Oil  is  used  at  -40*F. 

To  prevent  gas  leaks  from  occurring  on  cold  shots  utilizing 
the  threaded-motor  design,  motors  are  assembled  after  the  individual v: 
components  have  been  cooled.  The  clamp-type  motors  maybe  assembled 
at  ambient  temperature. 

3. 5  Instrumentation 
3.  5. 1  General 

The  basic  instrumentation  Bystem  used  for  both  .75C.50 
and  2C1.5-4  motors  is  shown  schematically  in  Fig.  21.  The  thrust  channel 
consists  of  a  dual^thrust  gauge,  filter,  calibrator,  amplifier,  and  an  analog 
and  digital  recording  system.  The  pressure  channel  is  basically  the  same 
as  the  thrust  channel  except  that  a  filter  is  not  used. 


II — r  ! — L,J~ — D> — 

THRUST  CALIBRATOR  FILTER  AMPLIFIER 
GAUGE 


C>- 

PRESSURE  CALIBRATOR  AMPLIFIER 

GAUGE 


ANALOG  RECORDING 
SYSTEM 


D,GlkliPPORD,NG 


ANALOG  RECORDING 
SYSTEM 


DIGITAL  RECORDING 
SYSTEM 


Fig  21  Thrust  and  pressure  instrumentation  used  with  micro-motors. 


'Consolidated  Electrodynamics  Corporation,  Pasadena,  California. 
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The  amplifiers  used  are  Kintel  Model  1121  which  have 
the  desirable  characteristics  for  transient  measurements  of  a  wide 
frequency  response  and  a  wide  linear  output -voltage  range.  A  wide 
frequency  response  is  not  in  itself  necessary  but,  unless  the.  cutoff  frequency 
of  the  amplifier  is  considerably  above  the  cutoff  frequency  of  the  filter,  it 
is  ncces sary  that  the  amplifier  have  the  characteristics  of  a  true  filter.  For 
example, '  amplifiers  characterized  by  a  slewing  rate  must  be  used  with  caution. 

Probably  the  best  types  of  integrating  systems  for 
transient  data  are  those  based  upon  the  voltage -to-frequency  type 
converter  such  as  the  Dymec  Model  221  IB.2  These  converters  have 
excellent  frequency  response,  wide  linear  overrange  capacity,  and  a 
full-scale  frequency  which  is  high  enough  to  provide  adequate  resolution 
for  shots  having  a  duration  of  1C  milliseconds  or  longer. 

3.5.2  Special  Considerations  for  .75C.50  Motors 

In  order  to  minimize  the  weight  of  the  overall  assembly, 
a  Bpecial  pressure  gauge  (Figs.  1  and  2)  was  designed  for  use  with. .75C.30 
motors.  This  pressure  gauge  is  the  same  basic  design  as  the  Rohm  &  Haas 
diaphragm  gauge  used  by  this  Division  for  a  number  of  years;  the  sensitive 
element  of  the  gauge  is  a  diaphragm  to  which  the  strain-gauge  bridge  and 
the  temperature-compensating  resistor  are  cemented.  The  gauge  internal 
cavity  is  filled  with  Celvacene  light  vacuum  grease  to  protect  the  diaphragm, 
reduce  the  motor  free  volume,  and  reduce  motor  heat  losses.  The  resonant 
frequency  of  the  assembly  using  a  50-pound  load  cell  is  approximately 
130  cycles/sec. 

The  motor,  pressure  gauge,  and  thrust  gauge  constitute 
a  highly  under -damped  second -order  system  which  is  shocked  into 
appreciable  oscillation  by  the  motor  thrust.  This  oscillation  is  of  such 
magnitude  that  it  interferes  with  the  digital  integration  by  causing  the 
thrust  gauge  output  to  swing  below  zero.  To  eliminate  this  oscillation, 
a  50-c.ycle  low-pa3s  filter  (Fig.  22)  is  included  in  the  thrust-gauge  circuit. 

rCohn  Electronics  Inc.,  Kintel  Div.  ,  San  Diego,  California, 

2Dymcc  Div.,  Hewlett-Packard  Co.  ,  Palo  Alto,  California. 


Fig. 22  Low-poss  filter  to  eliminate  thrust  oscillation  during  .75C.50 
motor  firings  . 


The  filter  is  designed  according  to  the  Butterworth1  condition  except 
that  R  is  adjusted  for  optimum  transient  response. 

The  attenuation  of  the  filter  above  the  cutoff  frequency 
increases  at  a  rate  of  nearly  18  decibe  1b /octave,  but  the  oscillations 
are  not  sufficiently  attenuated  when  a  50-pound  thruBt  gauge  is  used. 

To  provide  more  attenuation  the  capacitor  C  was  added  to  form  a  trap  for 
the  resonant  frequency  of  the  system,  approximately  130  cycles.  With 
this  addition  the  recorded  traces  are  oscillation-free. 

Thu  total  impulse  of  a  motor  is  defined  as/I^t)dt.  In 
practice,  however,  it  is  not  F(t)  that  is  integrated  but  a  function,  using 
Laplace  notation 

L*1A(s)B(s)F(s) 


'Frnst  A.  Guillemin,  *Synthesis  of  Passive  Networks,*  John  Wiley  fa 
Sons,  New  York,  1957. 
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where  F(s)  -  LF(t),  A(s)  is  the  transform  of  the  electronics  of  the  system 
and  B(s)  is  the  transform  of  the  thrust  stand  and  thrust  gauge.  If  the 
measured  impulse  is  to  be  valid  f  L-1A(s)B(s)F(8)dt  must  equal  /  F(t)dt 
even  though  L  1A(s)B(b)F(s)  only  approximates  F(t),  due  to  the  transient 
nature  of  the  thrust,  the  relatively  low  frequency  of  the  thrust  stand, 
and  the  filter. 

The  method  used  to  study  this  problem  was  to  assume 
as  the  thrust  function  F(t)  the  pulse  p(t)  -  p(t-a),  that  is,  a  square 
pulse  of  unit  amplitude  and  a  duration,  which  has  a  very  wide  frequency 
spectrum  and  to  evaluate  J*  L  1A(s)B(s)F(s)  by  integration.  The  frequency 
response  of  the  amplifier  and  the  digital  system  is  very  high  (10,000  cycles), 
and  the  only  frequency-dependent  component  of  A(s)  that  matters  is  the  filter. 
Let  A(a)  *s  PG(a)  where  G(s)  is  a  general  low-paBS  filter 

A  flm+A  +  ...  +1 

m  m-i 

- - - -  n  >  m 

B  sn  +  B  an"‘  +  ...  +1 
Yk  n-j 


If  B(s),  the  transform  of  the  thrust  stand,  is  ignored, 

and  since 

L[u(t)  -  u(t-a)]  =  ~  -  ~  e~aS 

B  O 

then 

.  A  sm  +  A  s1"'1  +  ,,,  +1 

LA(s)F(s)  =  -  (l  -  e"aS)  — - - p 

s  „  n  „  n-i 

Bs  +  B  s  +;..+! 
n  n-i 

By  the  theory  of  partial  fractions  the  right  hand  side  of  Eq.  2.  may  be 
expanded  into 

as,  jl_  ^  A  Be  +  C  Es  +  D  x 

[s  s  +  w,  (s  +  Pi)4  +  (s+pzP  +  U3  + 


(1) 

(2) 


P(l  -  e 


•  •  • 
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The  first  term  of  the  above  expression 


e"aS)  is  LF(t).  Therefore  if 


(i 


e 


-as 


) 


A 

8  +  Wj 


Ba  +  C  Es  +  D _ 

(a  +Pi)z  +  w|  (s  +  p2)*  +  <4 


=  0 
(3) 


then 


fm  ■  L 


L_lA(s)F(s)dt 


(4) 


That  Eq.  3  does  equal  zero  may  be  shown  by  direct  integration  although 
this  is  not  necessary  since, by  the  Second  Shifting  Theorem, L  e  afl<t>(s)  s 
<j»(t-a)u(t-a),  and  each  term  of  the  inverse  transform  of  Eq.  3  will  contain 
a  minus  exponential,  bo  the  integral  is  obviously  zero. 

The  analysis  could  be  extended  to  a  general  jT(t)  by 
dividing  F(t)  into  a  succession  of  pulses  Fj[u(t-iU)  -  u(t-L.-a)]  ani  then 
letting  a  **  0  . 

Since  the  thrust  stand  is  a  second-order  system  with  the 
familiar  equation 

mx  +  px  +  kx  “  F(t)  (5) 

its  tratteform  B(s)  is 
B 

msz  +  pe  +  k 


which  is  mathematically  the  transform  of  a  second-degree  filter. 
Therefore,  the  effect  of  the  thrust  stand  may  be  included  simply  by  adding 
another  term  to  Eq.  2  and  the  result  is  the  same.  Constants  are 
automatically  taken  care  of  by  the  calibration  procedure. 
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These  results  moan  that  it  is  possible  to  integrate  very- 
transient  forcing  functions  using  the  same  instrumentation  that  is  used 
for  larger  motors. 

To  demonstrate  the  mathematical  analysis  experimentally, 
a  series  of  .75C. 50-1.5  micro-motors  wae  fired  at  various  thrust/mass 
ratios  using  a  100 -pound  load  cell  and  various  ddghees'of  filtering.  The  precision 
of  the  measurements  was  not  good,  partly  because  the  severe  oscillations 
exceeded  the  range  of  the  load  cell.  However,  /Fdt  seemed  to  be  in¬ 
dependent  of  the  thrust/maas  ratio  and  filtering  (Table  V). 


Table  V 


Summary  of  Filter  Check  Data 


No. 

Rounds 

Batch 

Thruat/Masa^ 

Ratio 

Filter 

a 

*  1000 

5 

A 

10 

50  cycle 

234.8 

5 

A 

0.5 

15  cycle 

233.2 

6 

B 

10 

50  cycle 

236.0 

5 

B 

10 

None 

235.3 

4 

13 

0.5 

15  cycle 

235.5 

^See  Appendix  A  for  definition  of  ballistic  parameters. 
Thrust/mass  ratios  of  0.5  were  achieved  by  adding  50  lbs 
mass  to  the  standard  motor;  the  standard  configuration 
had  a  thrust/mass  ratio  of  10. 


In  general,  the  accuracy  of  the  measurements  (as 
distinct  from  the  problem  of  motor  reproducibility)  should  be  essentially 
the  same  as  with  larger  motors.  An  advantage  is  the  ease  with  which 
low -range  gauges  can  be  calibrated  with  dead  weights.  Occasionally, 
some  of  the  50-  and  ! 00-pound  thrust  gauges  have  given  trotible  because 
of  calibration  instability.  It  is  necessary  to  load  the  thrust  gauges  very 
conservatively  Sj.ii.ce  the  high  oscillation  level  can.  over-range  the  gauge. 


:5 


Sf 
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Generally  the  capacity  of  the  thrust  gauge  should  be  at  least  twice  the 
maximum  thrust  expected. 

Data  from  three  ten-motor  groups  of  .75G.50-3.5 
micro-motors  containing  RH-P-112cb  propellant  are  shown  in  Table  VI 
to  illustrate  the  precision  of  specific-impulse  measurement  which  can 
be  obtained  when  proper  experimental  techniques  are  used. 


Table  VI 


Specific- Impulse  Data  from  ,75 G .50-3.5  Micro-Motors 


Group  No. 

Group  Size 

Average 

FlMO 

1  <r 

(absolute) 

1  <r 

(per  cent) 

1 

10 

241.0 

1.0 

0.4 

2 

10 

241.3 

0.9 

0.4 

3 

10 

240.0 

0.4 

0.2 

3.  6  Data  Reduction 


Except  for  somewhat  different  definitions  of  parameters,  data- 
reduction  methods  for  micro-motors  are  the  same  as  those  for  large 
motors.  Conventional  data-r eduction  procedures  used  by  this  Division, 
as  well  as  the  exceptions  required  for  micro-motors,  are  given  in 
Appendix  A. 

Because  ignition  of  micro-motors  is  so  critical,  it  is 
possible  to  obtain  forms  of  pressure -time  traces  which  deviate  considerably 
from  those  of  larger  motors.  Application  of  conventional  data-r eduction 
techniques  to  pressure  traces  can  result  in  large  errors.  Several 
exaggerated  forms  of  pres  sure -time  traces  which  might  be  obtained 
with  micro-motors  are  given  in  Appendix  A. 

The  lower  trace  in  Fig.  23  provides  an  example  of  how 
erroneous  data  might  be  obtained  using  conventional  data-r  eduction  techniques.  The 
conventional  definition  of  start  of  action  time  and  burning  time  is  the  time 
at  which  pressure  has  risen  to  100  psia  on  the  primary  rise  portion  of 
the  trace  (point  3).  This  is  not  valid  for  the  trace  shown,  since  the  pressure 
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Fig  2.3  Desired  "square  trace  "(top)  and  an  unusual  pressure  trace  (bottom). 


up  to  point  4  is  caused  by  the  igniter  alone.  Propellant  begins  to  burn 

at  point  4,  and  this  should  be  used  as  the  start  of  burning  time.  If 

beginning  of  action  time  is  defined  as  the  time  at  which  the  thrust  trace 

begins  to  rise,  this  should  correspond  to  point  5  (the  time  at  which  the 

closure  ruptures).  Use  of  the  100-psia  convention  would  yield  unrealistically 

long  values  for  burning  time  and  action  time,  and  correspondingly  low 

values  for  r^e  P  ,  and  P^.  Values  of  F)000  would  be  too  high,  since  the 

low  value  of  would  be  used  in  the  C  equation  employed  in  the 

Test 


correction. 
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Tlie  most  desirable  pressure  traces  are  regular  as 
shown  in  Fig.  23;  if  enough  rounds  are  available  from  individual 
propellant  batches  the  ignition  system  can  usually  be  adjusted  to  give 
good  traces.  Pdor  traces  are  then  discarded  if  enough,  good  traces  are 
available.  However,  when  only  a  few  motors  are  available,  all  may 
produce  irregular  pressure  traces.  In  this  case  the  data  would  be  re¬ 
duced  using  procedures  outlined  in  Appendix  A, 

Any  correction  for  excessive  igniter  contribution  should 
be  made  once  F?c0o  *3  determined  by  the  procedures  outlined  in  Appendix  A. 
Defining  a  new  term  {Ffooe)i  as  the  igniter -corrected  F?00©,  the  following 
equations  may  be  derived! 


(FioooJj  =  Ftao 

(F?ooo).  *  Fjooo 


w. 


I  .  X  — - 
spi  w 


(for  correction  to  zero  igniter  weight) 


(6) 


I  1-1  -  0.0051  {for  correction  to  %%  igniter  weight) 


Bpi  I W 


(7) 


where  I  -  specific  impulse  of  igniter  (I  .  «  110  for  M-l  squib  and 
“P1  RIP  or  RHim)  8pi 

w^  *  weight  of  combustible  material  of  igniter 

w  =  weight  of  propellant 
P 

Measurement  of  igniter  specific  impulse  by  firing  igniters 
in  motors  containing  dummy  grains  and  measuring  /Fdt1  is  not  practical 
for  micro-motors  since  the  /  Fdt  of  0,05  grams  of  igniter  material  is 
not  measurable  on  a  50-lb  load  cell.  Consequently,  igniter  impulse  was 
determined  by  firing  multiple  rounds  with  various  amounts  of  igniter. 

The  resulting  plot  of  Fjooq  versus  igniter  weight  (Fig.  24)  was  extrapolated 
to  zero  igniter  weight,  yielding  (Fj000).,  The  specific  iripwlse  of  the 
igniter  was  then  deduced  by  comparing  with  F^ClOC  for  a  particular 

igniter  weight. 


*JANAF  Solid  Propellant  Rocket  Static  Test  Panel,  'Recommended  Procedures 
tor  the  Measurement  of  Specific  Impulse  of  Solid  Propellants,  M  SP1A 
Publication  No.  SPSTP/lOB,  March  1962. 
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NOTES: 

I  EACH  POINT  IS  A  FIVE-ROUND  AVERAGE 
2.  BAR  HEIGHT  DENOTES  RANGE  OF  F°100O 
FOR  THE  FIVE  ROUNDS 


0.2 

IGNITER  WEIGHT  -  gm. 


Fig. 24  Specific  Impulse  as  a  function  of  igniter  weight. 


The  standard  micro-motor  grain  length  is  inch 
shorter  than  the  case.  Occasionally,  however,  shorter  grains  will  be 
encountered,  either  because  insufficient  propellant  was  available  to  fill 
the  motor  or  because  the  cured  grains  were  trimmed  to  remove  end 
flaws.  Since  short  grains  will  expose  more  of  the  motor  case,  heat 
losses  will  be  higher  and  specific  impulse  correspondingly  lower.  It 
is  generally  desirable  to  correct  specific  impulse  for  the  additional 
heat  loss,  giving  a  specific  impulse  which  would  have  been  measured 
with  a  standard  grain.  This  correction  may  be  accomplished  by  use 
of  the  following  equation,  which  was  derived  from  equations  given  in 
Section  4  of  this  report. 


(F?0oo)i  -  ~J( F?ooo)2  +  fiq  gn 


actual 


I  nominal 


8.~a»iteJt>iWfanwteig  -<ttrft~aai»ar  tra«?:Mier  j  jaiiaBw 
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where  {F?000)jf 
J 

q 

«0 

n 

t 

c 

A 

m 


fJooo  which  should  be  obtained  with  standard  length  grain 
energy  conversion  factor,  =  778.16  ft-lbf/Btu 
heat  flux,  Btu/eec-ft2 

masB  conversion  factor  =  32.174  ft-lbm/lbf-sec2 
theoretical  thermal  efficiency 

time  interval  during  which  heat  is  lost  from  propellant 
gases  to  motor  hardware,  sec 

exposed  motor  hardware  surface  area,  in2 

propellant  weight,  Ibm 


Data-reduction  practices  used  to  obtain  P-K-r  data  from 

micro-motors  are  the  same  as  those  used  throughout  the  industry,  if 

a  regular  trace  is  obtained.  For  the  more  peculiar  traces,  the  procedures 

outlined  in  Appendix  A  for  micro-motors  Bhould  be  used  to  obtain 

P,  ■  r,  and  P  -  K  data, 
b  b  eq 

Due  to  the  small  web  thickness  (%  Inch)  of  .75C.50  motors, 
small  dimensional  variations  along  the  grain  havo  a  pronounced  effect 
on  the  apparent  burning  rate.  Any  variations  present  usually  cause  the 
propellant  to  burn  out  unevenly,  resulting  in  lower -than-normal  integral 
ratios  {/ P^dt// P^^^dt  as  defined  in  Appendix  A).  In  general,  the  lower 
the  integral  ratio  the  higher  the  apparent  burning  rate.  Therefore,  for 
burning-rate  measurements,  no  trace  is  considered  whose  integral  ratio 
is  not  at  least  0.90,  and  for  the  most  accurate  rates,  integral  ratios  of 
0.95  or  higher  are  preferred. 

Data  reduction  for  temperature-coefficient  firings  is 
handled  in  the  same  way  a3  for  P-K-r  data.  The  propellant  temporature 
upon  firing  is  estimated  from  experimental  curves  such  aa  those  in 
Fig.  21. 


I 

l 


Temperature  coefficient  of  pressure  at  constant  K,  t^., 
is  calculated  as  follows: 


whore  Pj  =  average  pressure  over  burning  time  at  temperature  Tj. 

P2  =  average  pressure  over  burning  time  at  temperature  T2. 

4  Data  Interpretation 

4.  1  Specific-Impulse  Data 

Strict  compliance  with  all  the  methods  and  techniques 
described  in  Section  2  of  this  report  will  yield  numbers  which  can  properly 
be  called  F°1000  in  the  (.75C.50-1.5,  .75C.50-2.5,  .75C. 50-3.5,  2C1.5-4) 
micro-motor.  These  numbers  in  themselves  have  no  particular  significance 
unless  they  form  the  basis  for  a  predicted  specific  impulse  for  the  specific 
composition  ip  a  much  larger  motor  which  has  a  specific  application, 

4.1.1  Theoretical  Scaling  Considerations  and  Equations 

Theoretically,  the  corrected  specific  impulse  of  a  motor - 
propellant  system  is  dependent  upon  the  following  factors: 

(a)  the  basic  available  combustion  energy  of  the  propellant; 

(b)  the  amount  of  energy  lost  from  the  propellant  gases  in  the  form  of 
heat  transfer  to  the  motor  components; 

(c)  the  degree  of  kinetic  and  thermal  equilibrium  attained  between  the 
gaseous  and  condensed  phases  of  the  combustion  products  during  the 
expansion  process; 

(d)  combustion  efficiency,  or  the  relative  completeness  of  combustion 
achieved  in  the  motor  chamber;  and 

(e)  the  degree  of  dissociation  and/or  recombination  of  the  product 
gases  during  the  expansion  process. 
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Items  (b)  through  (e)  can.  be  expected  to  vary  with  motor  size  and  con¬ 
figuration,  and  these  variations  are  responsible  for  the  frequently 
observed  specific -impulse  scaling  effect,  or  the  difference  in  F5mo 
between  motors  of  different  sizes  containing  identical  propellants. 

4. 1.1.1  Heat-Loss  Effects 

Assuming  complete  kinetic  and  thermal  equilibrium 
between  gaseous  and  condensed  combustion  products,  the  specific  impulse 
equation  may  be  written  as 

F!ooc  (AHci  =“\/~  (Hc  -He)  (10) 

where 

J  =  energy  conversion  factor,  778.16  ft -lbf/Btu 

gQ  ■  mass  conversion  factor,  32,174  ft-lbm/lbf-aec4 

AH  s  isentropic  enthalpy  change  of  combustion  pruducts  during 
expansion  through  the  nozzle,  Btu/lbm 

H  a  enthalpy  of  combustion  products  at  nozzle  entrance,  Btu/lbm 

.! '  ®  enthalpy  of  combustion  products  at  nozzle  exit,  Btu/lbm 

From  the  relation 

AH  a  CpAT  ,  (11) 


1  0  may  be.  rewritten  as 


where 


T^O 

*  1000 


1  - 


C  T  \ 
Pe  e 

"C  T" 

Pc  cl 


v  gn 


H  T) 
c 


(12) 


H 


specific  heat  at  constant  pressure  of  the  combustion  products 
(in  the  chamber)  at  temperature  Tc,  Btu/lbm- *R 

specific  heat  at  constant  pressure  of  the  combustion  products 
(at  the  nozzle  exit)  at  temperature  T  ,  Btu/lbm~*R 

chamber  temperature,  *R 


exit  temperature, 

C  'r 
Pe  c 


CPc  1  c 


C  T 
1\:  '■ 


•R 
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Assuming  that  heat  is  lost,  from  the  combustion  gases  in 

the  form  of  heat  transfer  to  the  combustion  chamber,  and  that  the  chamber 

enthalpy  would  have  been  H  instead  of  H  had  the  heat  loss  not  occurred, 

o  c 

the  heat  loss  may  be  expressed  as 


Q/m  =  H  -H  -  C  T  C  T  (13) 

o  c  p  O  -  P  c 

ro  rc 

where  Q/m  =  heat  lost  from  combustion  products  to  the  motor  chamber* 
Btu/lbm 

H  =  enthalpy  of  combustion  products  in  the  chamber  for  the 
case  of  no  heat  loss,  Btu/lbm 


specific  heat  at  constant  pressure  of  the  combustion 
products  at  temperature  Tq,  Btu/lbm-‘R 

chamber  temperature  for  the  case  of  no  heat  loss,  4R. 


Substitution  of  £q.  13  into  Eq.  12  yields 


^1000 


(H  -  U/m)q 
»-* 


-  Q/m 
^0 


(14) 


The  specific  impulse  that  would  have  been  obtained  for 
the  case  of  no  heat  loss  may  be  expressed  as 


C'  T» 

pc  e 

C  T 
p  o  • 


(15) 


where 


specific  heat  at  constant  pressure  of  the  exit  products  at 
temperature  T'  ,  Btu/lbm -*R 


AliJ  -  isentropic  enthalpy  change  of  the  combustion  products  during 

expansion  through  the  nozzle  for  the  case  of  no  heat  loss,  Btu/lbm 

H'  -  enthalpy  of  combustion  products  at  nozzle  exit  for  the  case  of 
no  heat  loss,  Btu/lbm 


■f 

e 

’1  1 


=  exit  temperature  tor  the  case  of  no  heat  loss,  *K. 
C  T' 

1  -  -LV-L 

C.  T 
!\>  " 
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Siin.ce  ^  a  r)lt  Eq.  15  may  be  substituted  into  Eq.  14 


to  yield 

F?ooo 


V 


l£  \z  .  2J  Q 
g0  ^  m 


Eq,  16  relates  the  specific  impulse  actually  obtained 
when  heat  loss  occurred  (Fjo^),  the  specific  impulse  which  would  have 
been  obtained  had  heat  loss  not  occurred  (-1®  ),  and  the  amount  of  heat 

rii-i 

lost  (Q/m). 

4, 1.1. 2  Two-Phaae-Flow  Effects 

The  specific  impulse  equation  may  be  written  as 


F?ooo 


mV+mV  mV  mV  yV+xV 

s  K  g  s  s  a  S 


b.,s  ■■  «  .--s-s ;  T  .  +  — ,  -s- 

pr+^rr  g0(mg  +  ™8)  «o  (mg  +  ™s]  go 


for  a  two -phase  system,  where 

m.  =■  mass  flow  rate  of  gaseous  phase  at  nozzle  exit,  lbm/sec 

S 

m  =  mass  flow  rate  of  condensed  phase  at  nozzle,  exit,  lbm/sue 
8 

V  =  velocity  of  gaseous  phase  at  nozzle  exit,  ft/aec 
g  . 

V  e  velocity  of  condensed  phase  at  nozzle  exit,  ft/aec 


(16) 


(17) 


x  ~ 


m  +  m 


rn 

i-Jr- 
m  +  m 


If  V^  =  (no  velocity  lag  of  the  condensed  phase). 


Eq.  I  7  reduces  to 


V 


VL 


(18) 
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Dividing  Eq.  17  by  Eq.  18  results  in 


o  yV  +  xV  V 

Llpm  =  — £ - 1  =  v  +  x 

V  y  *  V 

vL  g  g 

The  fractional  velocity  lag  may  be  defined  as 


V  -  V 


Substitution  of  Eq.  20  into  Eq.  19  yields 


— s  y  +  x(l  -  A)  =  y  +  x  -  xA  =  1  -  xA 

•*-1  «*  t 


Eq.  21  relates  the  specific  impulse  actually  obtained 
when  velocity  lag  of  the  condensed  phase  occurs  (F?^),  the  specific  impulse 
which  would  have  been  obtained  had  no  velocity  lag  occurred  (1®^).  the  W0l8ht 

fraction  of  condensed  products  (x),  and  the  fractional  velocity  lag  of  the  con¬ 
densed  products  (A). 

4.  1.  1.3  Combined  Heat  Loss  and  Two-Phase  Flow  Effects 


*HL 

ptr 


Eq.  21  maybe  written  in  the  form 


1  -  xA 


where  =  specific  impulse  which  should  be  obtained  if  velocity  lag 

occurs  but  heat  loss  dpes  not  occur 

I00  =  specific  impulse  which  should  be  obtained  if  neither  heat 

loss  nor  velocity  lag  occurs. 

Substitution  of  Eq.  22  into  Eq.  16  yields 


(I00)2  (1  -  xA)2  - 


2  J  Q 


i 

II 


(f10»)2  =~\f  f  I  x'AJ|  ^Fiooo)i  +  (Q/m)!  -  |£f|(G/xn)j 


The  heat-loss  term,  Q,  in  Eq,  26  may  be  expressed  as 


Q  * 


qt  A 
z 

144 


where  q  =  heat  flux,  Btu/sec-ft2 

t  =  time  interval  during  which  heat  is  lost  to  internal  motor 
hardware  (motor  action  time),  sec 

A  *  surface  area  of  internal  motor  hardware  to  which  heat 
ia  last,  in* 

Substitution  of  Eq,  27  into  Eq.  26  yields 

2/1 ms 

V  |1  -  xAt 


(F?ooo>* 


(F?ooo)?  +  144^ 


t  AM 

2  Jr|  q  1 

t  A 

z 

2) 

‘  144  g()  ! 

m 

2 
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The  unknown  quantity  I00  may  be  eliminated  by  writing  Eq.  23  for  two 
static  test  motors  and  solving  the  resulting  equations  simultaneously  as 
follows : 

(Ffoootf  =  (I00)2  (1  -  xAtf  -  ZLV  (Q/m)j 


(F?ooo)l  =  (I00)2  (1  -  xAj)2  -  (Q/m)2 


(24) 

(25) 


(27) 


(28) 


Eq.  28  relates  the  specific  impulses  of  two  arbitrary 
static  test  motors  (containing  identical  propellants)  in  terms  of  their 
respective  velocity  lags  and  heat  losses,  and  is  the  basic  equation  used 
by  this  Division  to  predict  large-motor  specific  impulse  and  to  determine 
scaling  factors  from  small-motor  static-test  data. 


s 

m 


1 


I 
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Eq„  28  is  not  a  rigorous  equation;  many  simplifying 
assumptions  were  made  in  its  derivation.  However,  it  has  proven  to 
be  valuable  and  useful  in  its  present  form,  and  consequently  no  attempts 
have  been  made  to  perform  a  more  rigorous  analysis. 

4. 1 . 2  Calculation  of  Fractional  Velocity  Lag  (A) 

A  FORTRAN  computer  program  developed  at  this  Division 
is  currently  used  to  calculate  fractional  velocity  lag  (A).  A  simplified 
mathematical  model  is  used  to  approximate  the  real  flow  process.  This 
simplified  model,  which  ignores  the  coupled  effect  of  the  particle  lag 
on  the  flow  of  the  gas  phase,  permits  rapid  estimation  of  the  effects  of 
solid  particles  and  nozzle  configuration  on  performance.  The  model  is 
not  sophisticated;  however,  it  has  been  satisfactory  for  the  determination 
of  specific -impulse  scaling  factors  from  small-motor  static  tests. 
Description  of  the  program  and.  the  mathematical  model  is  given  in 
Appendix  £5  and  a  complete  listing  of  the  FORTRAN  program  ia  given 
in  Appendix  C. 

Before  the  present  program  was  developed,  an  even  more 
simplified  version  programed  for  a  Royal-McBee  LGP-30  computer  was 
used  to  calculate  values  of  A,  For  a  given  case,  the  current  program 
sometimes  yields  A  values  quite  different  from  those  of  the  original 
program.  Consequently,  scaling  factors  derived  from  a  given  set  of 
small-motor  data  could  differ  somewhat,  depending  on  the  program  used 
to  calculate  A  values.  However,  practically  identical  predictions  of 
specific  impulse  in  a  given  large  motor  are  obtained  with  both  velocity- 
lag  programs. 

Use  of  A  values  calculated  for  our  reference  RH-P-112 
plastisol  propellant  generally  result  in  valid  predictions  for  other  propellant 
compositions,  even  though  the  derived  scaling  factors  might  differ 
considerably  from  those  obtained  if  A  values  calculated  for  the  specific 
propellant  are  used.  F'or  this  reason,  a  set  of  A  curves  generated  for 


$ 
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RH-P-112  propellant  are  shown  in  Appendix  D.  Curves  of  A  vs  nozzle 
diameter  at  various  nozzle  expansion  ratios  and  particle  diameters  are 
shown  for  the  Rohm  &  Haas  .75C.50,  2C1.5,  and  6C5  motor  nozzles,, 
for  the  USAF  BATES  motor  nozzle,  and  for  the  Aerojet  10KS-2500 
motor  nozzle.  These  curves  will  generally  result  in  valid  specific 
impulse  predictions  for  most  propellant^  but,  if  a  computer  is  available, 
use  of  A  values  calculated  for  the  actual  propellant  is  recommended. 
Additional  computer -calculated  A  values  will  obviously  be  required 
for  predictions  in  motors  other  than  those  for  which  curves  are  given 
in  Appendix  D. 

4.  1 . 3  Determination  of  Scaling  Factors 

Examination  of  Eq.  28  indicates  that  an  (Ffooo)2  for  some 

arbitrary  motor  could  be  calculated,  baBed  on  the  ballistic  results 

obtained  from  static  teats  of  some  other  motors,  if  values  for  A  and  q 

■were  known.  The  weight  fraction  of  condensables,  x,  is  known  from 

thermochemical  calculations;  values  for  and  m,  are  known 

z 

from  static  testa;  A,  and  Aa  may  be  calculated  on  the  basis  of  the  respective 
motor  configurations;  and be  closely  estimated,  based 
on  the  ballistic  results  from  motor  1  and  the  configuration  of  motor  2, 

Further,  by  the  procedure  described  in  Appendices  B 
and  C,  only  the  particle  diameter  (N)  must  bo  specified  to  calculate 
Ai  and  A2.  All  other  parameter  j  necessary  for  the  fractional  velocity 
lag  calculation  are  known,  once  the  propellant  composition  and  nozzle 
configurations  have  been  specified.  This  means  that  (F?000)2  could  be 
calculated,  based  on  static  test  data  from  motor  i,  if  N  and  q  were  Jcnown. 

Many  simplifying  assumptions  were  made  in  deriving 
Eq.  28.  Probably  the  most  questionable  is  the  assumption  that  specific 
impulse  scaling  is  due  only  to  boat  loss  and  two- phase-flow  effects. 
Consequently,  it  is  doubtful  that  use  of  the  real  values  of  N  and  q  in 
E.q.  28,  even  if  these  values  could  be  directly  measured,  would  provide 
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proper  correlation  of  the  specific  impulses  of  different  motor  configurations. 
For  this  reason,  effort  was  directed  toward  a  determination  of  the 
effective  or  apparent  values  of  N  and  q  which  embody  all  questionable 
assumptions  but,  although  they  may  differ  considerably  from  the  real 
values,  nevertheless  provide  correlation  between  specific  impulses 
of  different  motor  configurations.  These  effective  values  of  N  and  q 
are  called  scaling  factors. 

By  assuming  that  the  heat  flux  (q)  is  independent  of 
motor  configuration,  Eq.  28  may  be  solved  for  q  to  yield 


The  utility  of  Eq.  29  i*  immediately  obvious.  If  static 
test  data  are  available  in  two  motor  configu.  l-.ons,  Eq.  29  is  an  equation 
in  the  two  unknowns  q  and  N,  If  static  test  data  are  available  in  three 
motor  configurations,  a  total  of  three  aquations  in  the  form  of  Eq.  29 
may  be  written  in  the  unknowns  q  and  N;  this  is  more  than  sufficient  to 
allow  simultaneous  solution  for  q  and  N. 

Since  the  value  N  does  not  appear  explicitly  in  Eq.  29$ 
the  method  of  intersecting  curves  must  be  used  to  effect  a  solution.  All 
quantities  derived  from  static  teats  and  thermochemical  calculations  are 
substituted  in  Eq.  29  for  the  two  motors  to  be  compared.  Various 
values  of  N  are  assumed  and  values  of  A.,  and.  q  are  calculated  for 
each  until  an  extensive  curve  of  q  versus  N  has  been  generated.  A 
different  combination,  of  two  motors  is  selected  and  another  curve  of  q 
versus  N  is  generated.  These  curves  will  intersect  at  the  point  of  solution 
for  the  scaling  factors  q  and  N,  and  theso  scaling  factors  may  then  be 
used  in  conjunction  with  Eq.  28  to  predict  specific,  impulse  in  any  other 
motor  configuration. 
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Static  teat  data  must  be  available  from  a  minimum  of  three 
motor  configurations  to  allow  solution  for  q  and  N  by  the  above  method, 
but  more  accurate  scaling  factors  will  result  if  data  are  available 
from  more  than  three  configurations.  Also,  more  reliable  scaling  factors 
should  be  obtained  from  large  motors  than  from  small  motors,  due 
primarily  to  the  difficulties  associated  with  accurate  measurement  of 
ballistic  parameters  in  small  motors.  Full  use  should  be  made  of  the 
maximum  quantity  of  propellant  available. 

For  best  results,  the  motors  used  to  obtain  the  data 
necessary  for  scaling-factor  determination  should  be  of  two  distinct 
types:  (1)  motors  having  the  same  grain  O.  D.  and  grain  cross -Beetional 
configuration  but  differing  in  length  by  at  least  a  factor  of  two,  and  (2) 
motors  which  are  linearly  scaled  but  which  differ  in  grain  O.  D.  by  at 
least  a  factor  of  two.  The  constant -O.  D. ,  different -length  motors  will 
have  similar  two-phase-flow  losses  but  will  differ  significantly  in  heat 
losses,  while  the  scaled  motors  will  have  practically  the  same  heat  loss 
per  unit  mass  (Q/m)  but  will  differ  -significantly  in  two -phase -flow  losses. 
Thus,  any  F^ooo  differences  for  the  different-length  motors  should  be 
primarily  dependent  on  heat  losses,'  whereas  FJomj  differences  for  scaled 
motors  should  be  primarily  dependent  on  two-phase-flow  losses.  The 
q  versus  N  curve  for  the  scaled  motors  should  then  be  much  steeper  than 
the  q  versus  N  curve  for  the  different-length  motors,  providing  a  very 
definite  intersection  point  for  scaling  factor  determination 

The  micro-motor  system  described  in  this  report  has 
been  successfully  used  to  determine  scaling  factors  for  propellants 
available  in  very  limited  quantities.  The  minimum  three -motor  con¬ 
figurations  used  are  the  .75C.50-1.5,  the  .75C.50-3.5,  and  the  2C1.5-4. 

Due  to  the  difficulties  associated  with  measurement  of  ballistic  parameters 
in  such  small  motors,  multiple  firings  are  generally  conducted,  using  a 
minimum  of  five  .75C.50-1 ,5,  five  .75C.50-3.5,  and  three  2C1.5-4  motors 
if  possible.  These  tliirteen  rounds  can  be  cast  from  a  nominal  2 -lb 
propellant  mix. 
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Scaling  factors  can  be  estimated,  based  on  a  knowledge 
of  scaling  factors  for  similar  propellants,  if  they  cannot  be  determined 
experimentally  because  of  propellant  limitations.  Estimated  scaling 
factors  will  be  questionable  at  best,  but  may  be  of  some  use  in  making 
approximate  scale-up  predictions. 

If  static  teat  data  are  available  only  for  one  motor 
configuration,  both  q  and  N  must  be  estimated.  The  following  comments 
are  offered  as  a  guide  in  making  these  estimations. 

Results  from  a  series  of  tests  of  plastiaolnitrocelluloae  propellant 
having  aluminum  concentrations  from  5  to  22%  by  weight  were  correlated 

r 

using  Eq.  29  to  determine  scaling  factors.  For  all  compositions  tested, 
q  values  fell  between  1100-150Q  Btu/sec-ft*  and  N  values  between  1,1- 
1.8  microns.  Several  NFPA  compositions  having  aluminum  concentrations 
from  10-15%  gave  q  values  in  the  same  range  and  N  values  from  0.8- 
1.1  microns.  The  q  values  could  be  roughly  correlated  by  tho  equation 

q  =  6.382  (T  /1000)4  +  156.4  (T  /1000)  (30) 

c  c 

whore  T  =  theoretical  flame  temperature,  *K.  The  empirical  constants 
c 

in  Eq.  30  were  determined  using  a  heat  flux,  q,  of  1300  Btu/sec-ft*  for 

RH-P-112cb  propellant  (T^  =  3325*K),  ai*d  assuming  that  the  heat  flux 

was  60%  radiative  (proportional  to  T^4)  and  40%  convective  (proportional 

to  T  ). 
c' 

If  data  are  available  for  only  one  motor  configuration, 
an  estimated  q  should  be  calculated  from  Eq.  30.  An  estimated  value 
for  N  will  have  to  be  arbitrarily  selected,  based  on  considered  judgment. 

If  data  are  available  for  two  motor  configurations  only,, 
a  curve  of  q  versus  N  should  be  constructed  as  previously  outlined.  A 
value  for  q  should  be  calculated  from  Eq.  30,  and  the  N  value  which 
corresponds  to  the  calculated  q  should  be  picked  from  the  q  versus  N 


curve. 
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A  detailed  example  illustrating  scaling  factor  determination 
from  .75C.50-1.5,  .75C  50-3.5,  and  2C1.5-4  motor  static  test  data  is 
presented  in  Appendix  E.  Auxiliary  curves  permitting  rapid  determination 
of  the  area  term  (A)  for  selected  motors  in  Eqs.  28  and  29  are  given 
in  Appendix  F. 

4. 1 . 4  Prediction  of  Large -Motor  Specific  Impulse 

Large -motor  specific  impulse  predictions  are  made 
using  Eq.  28  and  are  relatively  straightforward  once  the  scaling  factors 
have  been  determined.,  The  configuration  of  the  large  motor  is  selected, 
allowing  determination  of  t2,  A2,  m2,  and  A2.  These  parameters  are 
substituted  into  Eq.  28  along  with  the  corresponding  parameters  for  one 
of  the  small  motors  for  which  static  test  data  are  available,  and  (F^o^o^ 
is  calculated  for  the  large  motor. 

A  detailed  example  illustrating  the  prediction  procedure 
is  given  in  Appendix  E. 

Micro-motors  should  not  be  used  to  make  direct  specific 
impulse  comparisons.  That  is,  specific  impulses  of  a  test  propellant 
and  a  reference  propellant  as  determined  in  micro-motors  Should  not  be  ex¬ 
trapolated  to  some  large  motor;  biased  results  will  almost  certainly 
be  obtained  unless  the  propellants  being  compared  have  similar  heat 
losses  and  two-phaserflow  properties.  Flame  temperatures,  burning 
rates,  and  concentrations  of  condensable  products  in  the  exhaust  must  be 
about  equal  for  the  propellants  before  comparative  specific  impulses 
have  any  significance. 

4.  2  P-K-r  and  Temperature-Coefficient  Data 

While  specific-impulse  values  obtained  in  micro-motors 
must  be  scaled  to  larger  motors, P-r  and  tt  data  obtained  in  micro-motors 

XY. 

are  t.he  same  as  for  large  motors,  provided  that  proper  experimental, 
techniques  are  used  in  making  the  measurements.  A  scaling  effect 
exists  for  P-K  dntn,  however. 
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The  P-K-r  relationships  between  micro-motors  and 
larger  motors  were  evaluated  by  firing  a  plastisol  nitrocellulose  com¬ 
position  in  fourteen  different  motor  configurations.  The  motors  ranged 
in  size  from  a  .75C.50-1.5  micro-motor  (10  gm)  to  a  9C7.5-17  (19  lb)  motor. 
Data  from  three  complete  series  are  tabulated  in  Table  VII.  The 
propellant  burning  rate  was  varied  by  vising  ammonium  perchlorate 
having  different  particle  sizes. 

The  values  of  K  in  Table  VII  were  obtained  by  correcting 
the  actual  values  to  1000  psi  (assuming  n  =  0.5)  so  that  a  comparison 
could  be  made  between  motor  sizes.  This  correction  was  small  in  most 
cases,  since  all  firings  had  been  designed  to  give  an  average  pressure 
of  1000  psi.  At  mass-discharge  rates  greater  than  3  lb/sec  K  was 
essentially  constant  but  at  mass-discharge  rates  of  less  than  3  lb/sec 
K  decreased  with  increased  mass  discharge  (Fig.  25).  The  probable 


Table  VII 

Ballistic  Data  from  F '  1 •  of  Vgrlom-31gQ  Motors 
Containing  Plaptltol  Propellants  Having  Thrco  Different  Burning  Raton 


Grind  "n»  APCa _  Grind  *b»  APCft _  Qrlr.d  APC* 


Motor 

Designation 

K 

TTI 

Tl> 

iVi 

0.094 

K 

trt 

rb 

(in/sec) 

USZliSl 

K 

. 

rh 

(!**/•  ■«) 

ih 

(lb/sec) 

.75C.50-1  .5 

1  78.0 

0.70 

213.9 

0.57 

0,06b 

264,4 

0.43 

0.061 

.75G.50-I.5 

177.5 

0.70 

0,15 

--- 

— 

.... 

.750,50-3.5 

I  .’6.8 

0.67 

0.22 

---- 

*  ■ 

•*  - 

,750.30-9 

171.6 

0.61 

.... 

— 

---- 

*  • 

— 

-SCI. 5-2 

175.4 

0.71 

0.29 

21  1.4 

0.56 

0.25 

273.1 

0.47 

0.21 

iC  1.5-4 

173.0 

0.6ft 

0.67 

?i  1.4 

0.56 

0.53 

263.4 

0.45 

0.45 

2  Cl  .5  -7.5 

170.7 

0.67 

1.29 

205.9 

0.56 

1.02 

265.1 

0.45 

0.04 

2CI.5-II  ,i 

166.3 

0.70 

1.99 

206,8 

0.57 

1.59 

259.4 

0.45 

1.30 

2C1.5-I  -l.b 

166.2 

0.69 

2.6  1 

205.1 

0.56 

2.15 

253.8 

0.4-7 

'.73 

6C*1  - 1» 

165.3 

0.69 

3,60 

198,5 

0.56 

2.92 

24  3,1 

0.46 

2.15 

6C5-I  l.-l 

164.6 

O.bfl 

7.0  b 

196.2 

0.55 

5.75 

243.1 

0.45 

■1.68 

hCr--24 

164.5 

0.7* 

1  1.38 

197.1 

0.55 

1  l.bR 

248.1 

0.45 

9.71 

6C5-3  3 

161.3 

0.69 

1  9  .Jto 

1  99,4 

0.55 

1  5.P.0 

247.7 

0.45 

!  3. -12 

•K’.  i.S-l  7.  1 

U.l.'l 

!  r«.*<ii 

l-m.-i 

n.'ih 

12.60 

247.7 

0.45 

1 0.59 
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Fig.  25  Kjooo  05  0  ^unc^on  ro  ar*d  ammonium  perchlorate  particle  size 
for  plaslisol  proptfliuni . 

cause  Is  a  combination,  of  heat  loss  and  poor  combustion,  efficiently.  Curves 

of  K  af-  r  ction  of 
rn 


a  heatrlost  in  which  t  is  burning  time,  A  is  area  of  hardware 

exposed,  a  ;e  propellant  weight  are  shown  in  Fig,  26  for  propellant 

containing  :  a"  ammonium  perchlorate.  A  different  curve  was 

obtained  fv  an  motor  diameter,  indicating  that  heat-loss  variations 
alone  are  not  responsible  for  the  change  in  .  As  motor  length  in¬ 
creased  K  tended  to  become  constant. 


1 


I 


II  Hill  i 


The  burning  rates,  for  the  fourteen  motor  sizes  in 
each  series  were  corrected  to  1000  psi  (Table  VII).  These  burning 
rates  were,  within  experimental  error,  independent  of  mass  discharge 
rate  or  propellant  web  thickness. 

The  results  indicate  that  micro-motor  P-K-r  data 
can  be  very  useful  on  new  propellant  compositions.  Burning-rate  data 
can  be  used  directly,  and  P-K  data  can  be  used  after  the  proper  cor¬ 
rections  have  been  made. 

A  pressure  value  corresponding  closely  to  the  equilibrium 
pressure  of  the  rocket  motor  would  be  most  meaningful  in  P-K  relation¬ 
ships.  The  average  pressures,  P  for  two  motors  having  the  same 
amplitude  pressure  trace  but  different  burning  times  are  different  since 
the  aspect  ratio  of  the  pressure -time  trace  varies  with  motor  size.  Use 
of  points  5  and  7  on  the  pressure-time  trace  (see  Appendix  A)  as  the 
limits  of  the  equilibrium  pressure  integral  is  recommended.  This  pressure 
should  be  essentially  constant  for  a  given  K  regardless  of  the  initial  rise 
time  or  the  total  burning  time  of  the  propellant  and  should  be  ttsed  in 
P-K  calculations. 
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The  average  pressure,  P  ,  would  still  be  used  for  burning  rate  measurements, 
The  data  of  Tables  VELI  and  IX  illustrate  temperature 
coefficient  relationships  between  micro-motors  and  6C5-11.4  motors. 

Data  obtained  from  firing  a  series  of  .75C. 50-1.5  micro-motors  with  a 
typical  plastisol  propellant  at  the  extremes  of  ambient  temperature 
(35 'F  to  71  *F)  are  shown  in  Table  VIII.  Using  theBe  data,  it.. 

TV 

(35*F  to  71  *F)  =  0.21%/*F.  The  temperature  coefficient  for  this 
propellant  as  determined  in  the  6G5-11.4  motor  was  ir  (-40*F  to  135*F)  - 
0.23%/*F  while  the  temperature  coefficient  determined  in  insulated 
micro-motors  was  -rr  (-33*F  to  97*F)  =  0.25%/*F, 

Xv 


Table  VIII 

Temperature  Coefficient  Data  from  -75C.50-1.5  Micro-Motors 
at  Ambient  Temperatures 


Temperature 

TT 

,  *F  b 

K  s  164 

r .  fP.dt/fP^  t  ndt  P  /P. 

b  J  b  total  max  b 

37 

983 

0.601 

0.94 

1.19 

37 

997 

0.603 

0.95 

1.17 

37 

981 

0,601 

0.95 

1.19 

37 

1004 

0.590 

0.96 

1.16 

37 

983 

0,609 

0.95 

1.19 

37 

1005 

0.594 

0.96 

1.14 

Avg  37 

99  2 

0,600 

0.95 

1.17 

71 

1059 

0,643 

0.96 

1.13 

71 

1072 

0,654 

0.95 

1.15 

71 

1072 

0.672 

0.95 

1.13 

71 

1  068 

0.652 

0.94 

1.14 

71 

1063 

— 

0.92 

1.13 

71 

1  049 

---- 

0.92 

1.14 

Avg  71 

1  064 

0.65  6 

0.94 

1 .14 

,tk  - 

0.21%/T 
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Tablc  IX 

Tcmpc rat ure  Coefficient  Data  from  Insulated  .75C.5Q-1.5  Micro-Motors 


Temperature,  *F 
-33 
-33 
-33 
-33 

Avg 

97 

97 

Avg 


K  =  186 


^b 

rb 

758 

0.470 

717 

0.449 

820 

0.530 

779 

0.501 

769 

0.488 

1078 

0.705 

1061 

0.660 

1069 

0.682 

fp,  dt/r  „  at 
J  b  ' J  total 

0.93 

0.95 

0.92 

0.92 

0.94 

0.88 

0.92 

0.90 


irK  c  0.25%/ *F. 


These  results  indicate  that  good  temperature-coefficient 


data  can  be  obtained  in  micro -motors 
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Appendix  A 

Data  Reduction  Definitions  and  Procedures 

Typical  pressure-time  and  thrust-time  traces  are  shown 
in  Fig.  A-i.  While  actual  oscillograph  traces  have  heavy  vertical  bars 
every  100  msec,  lighter  vertical  bars  every  10  msec,  and  a  1000-cycle 
sine  wave  to  provide  1-msec  timing,  these  extraneous  features  have  been 
omitted  for  clarity.  Also  deleted  from  the  traces  is  a  train  of  "blips" 
which  indicate  cumulative  instantaneous  pressure-time  and  thrust-time 
integrals. 


Fig  A-I  Typical  pressu.e  lime  lux,'  (lop)  and  thrust -time 


trace  (bottom). 


A  -2 


The  following  time  points  in  Fig.  A-l  are  defined  as  follow*; 
tj  =  time  of  current  application  to  squib 
t2  =  time  of  initial  pressure  rise 
t3  =  time  when  pressure  has  increased  to  100  psia 
t4  =  time  when  propellant  begins  to  burn 
faj  »  time  when  thrust  beginc  to  rise 
t$  *  time  when  equilibrium  pressure  is  attained 
£7  B  time  of  end  of  burning 

tf  »  time  when  pressure  has  decayed  to  100  psia 
tj  -  time  when  pressure  has  decayed  to  zero  psig 
t|o  B  time  when  thrust  has  decayed  to  zero  lbf. 

The  following  definitions  and  formulas  describe  the 
calculations  performed  from  the  static-test  oscillograph  records  and 
from  the  digital  instrumentation  to  determine  ballistic  data. 

1.  Standard  Calculations  frcm  Oscillograph  Traces 

Zp(p«i/in)  *  pressure  factor  b 

(calibration  value  of  step)  X  (value  of  pressure  gauge,  pal) 
(height  of  step,  in)  X  (gaug'd  calibration  factor) 


t » 


I 

i 


Z_(lbf/in)  =  thrust  factor —obtained  in  the  same  way  as  Z_. 


de  "  Qiii|  (z4  -  zj 

s 

1 

(NOTE:  Numerical  subscripts  and  integral  limits  are  abbreviated  designations 
of  times  shown  in  Fig.  A-l.  Thus  F3  is  the  instantaneous  thrust  at  tj) 

Ignition  Delay(msec)  -■  t3  -  tj 

t  (msec)  =  equilibrium  pressure  time  *  tf-ti 
eq 

t  (sec)  =  burning  time  =  t7  -  t3 

_  /•  /  \  A  ,  web  (in)  XI 000 

r^in/sec)  =  average  burning  rate  over  burning  time  »  — — — (rngecy - 

D 
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t  (msec) 
ap 

laf  (msec^ 

t-  (in/sec) 
a  '  ' 


P 

P 

F 


,  (psia)  s 

t 

(psift) 

(Ibf) 


max 


max 


=  pressure  action  time  -  tg  - 

-  thrust  action  time  =  ti0  -  t8 

-  average  burning  rate  over  action  time 


web  (in)  X  1000 

**•  f  7 _  “  %  r 


t  ^  (msec) 

■  ignition  pressure  =  (in)  X  Zp  (pai/in)  +  14.41 

=  maximum  pressure  =  P  (in)  XZ_  (psi/in)  +  14.4 

max  r 

=  'maximum  thrust  =  F  (in)  X  Z_  (lbf/in) 

max  F 


lAverage  atmospheric^pressure  at  teat  facility. 


Z.  Special  Calculations  from  Oscillograph  Traces 

If  digital  data  are  not  obtained  or  are  unsatisfactory, 
the  areas  under  the  pressure -time  and  thrust-time  traces  are  manually 
integrated  with  a  pl&himeter  and  the  following  calculations  are  performed. 


Z  (msec/in)  time  factor 

t 

total  number  of  I-rasec  timing  marks  in  an  interval 
length  of  i.he  interval  "(in) 

F  (psia)  b  equilibrium  pressure 

Si  Pdt  (in2-)  X  Z  (msec/in)  X  Z  (psi/in) 

(tf  -  ts)  (msec) 

(psia)  =  average  pressure  over  burning  time 

/7Pdt  (in2)  X  Z  (msec/in)  XZp  (psi/in) 

*  -2 - L  -  -  +14  4 

(t7  -  t,)  (msec) 

P  (psia)  =  average  pressure  over  action  time 

a 

8 

/  Pdt  (in2)  X  Z  (msec/in)  X  (psi/in) 


=  _JL 


(t8  -  t3)  (msec) 


14.4 
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■^Ptotaldt  (Psia"8ec)  = 
P  Pdt  (In2)  X  Z 
=  J _ _ _ 1 


total  pressure-time  integral 

(msec/in)  X  (psi/in)  +  14.4  X  (t^  -  tj)  (msec) 


/Fdt  (lbf-eec)  =  thrust-time  integral 

/10Fdt  (in2)  X  (msec/in)  X  Z_  (lbf/in) 

s  J _ _ _ _ 

1000 


Fa  <lb£> 


*  average  thrust  over  action  time 


3 


Fdt  (lbf-sec) 
*io  -  tj)  (msec) 


X  1000 


3,  Calculations  from  Digital  Integrals 

3. 1  Calculation  of  Total  Integrals 

The  following  data  are  printed  out  from  the  Dyxnec  digital 
instrumentation  system  for  both  pressure  and  thrust  channels: 

Zero  (baseline),  counts/sec 
All  calibration  steps,  counts/sec 
Gross  integral,  counts 
Gate-open  time,  sec 

The  following  calculations  are  made  using  these  data: 

Net  full  scale  (c/sec)  *  gross  full  scale  (c/sec)  -  zero  (c/sec) 

Base  area  (counts)  *  gate  open  time  (counts)  X  zero  (c/sec) 

Net  integral  (counts)  »  gross  integral  -  base  area 
To]>Btep  pressure  value  (psi) 

(calibration  value  of  top  step)  X  (value  of  pressure  gauge,  psl) 
gauge  calibration  (actor 

Top-step  thrust  value  (lbf) 

(calibration  value  of  top  step)  X  (value  of  thrust  ,  gauge,  psi) 
gauge  calibration  factor 
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/  Fdt  (lb£-sec)  =  thrust -time  integral 

(top-step  thrust  value,  lbf)  X  (net  integral,  counts) 

( not  full  scale,  c/aec) 

J 'P  ^dt  (psia-sec)  =  total  pressure -time  integral 

(top-atep  pressure  value,  psia)  X  (net  integral,  counts)  , 
(net  full  scale,  c/aec) 

+  14.4  X  (gate  open  time,  sec) 


3.2  Calculation,  of  P  ,  P.  ,  and  P 
_ _ eq  b _ a 

P  (psia)  -  equilibrium  pressure  -  1a  +  14.4 

eq  (t7  -  t4)  (sec) 

_  /TPdt  (psi-sec) 

P  (psia)  =  average  pressure  over  burning  time  =  /?■  1  /  «— — -  +  14,4 

b  tt7  "  ls ;  l»ec; 

_  /^Pdt  (psl-sec) 

P&  (psia)  a  average  pressure  over  action  time  •*  — t  )  (neo)*  * 

The  respective  pressure-time  integrals  required  to 
calculate  P  ,  P^»  and  are  determined  in  a  manner  completely 
analogous  to  that  outlined  for  determining  the  total  pressure-time  integral. 
However,  the  respective  gross  integrals  in  counts  are  not  printed  out 
and  must  be  determined  by  counting  "blips'*  representing  the  instantaneous 
cumulative  integral.  These  blips,  generated  by  the  digital  instrumentation 
are  visually  displayed  on  the  oscillograph  traces  in  groups  of  5,  50,  and 
500.  Each  blip  is  equivalent  to  two  counts.  Also,  the  respective  proper 
time  intervals  in  sec  (t7  -  t$,  t7  -  t3,  tg  -  ts)  must  be  used  in  place  of 
gate -open  time. 


4. 


Subsequent'  Calculations 
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K.  =  S  /A  .  where  S  1  is  an  integral  average  surface 

area  and  A^,  is  the  arithmetic  average  of  throat  area 
before  and  after  firing.  Throat  diameter  after  firing 
is  measured  before  any  accumulated  deposit  is 
removed. 

XR  b  Nozzle  expansion  ratio  ■  (D^/D^)2,  where  and 

are  the  respective  arithmetic  average  values  of 

nozzle  throat  diameter  before  and  after  firing  and 

nozzle  exit  diameter  before  and  after  firing. 

C_  (lbm/lbf-sec)  =  Discharge  coefficient  =  c ^ e* 8 

ATX‘'Ptotal 

Charge  weight  is  the  difference  between  weight  before 
firing  {before  inhibitor  is  applied)  and  weight  after  firing. 
p  Discharge  coefficient  corrected  for  pressure  drop 
°  along  the  length  of  the  grain  =  Cjj(P^/Pc),  where  P^ 
is  the  head-end  pressure  and  is  the  stagnation 
pressure  at  the  entrance  to  the  converging  section  of 
the  nozzle. 

Calculate  Jq  (initial  ratio  of  throat  area  to  port  area) 
and  from  Fig.  A-2  find  the  corresponding  Mach  number  (M  ).  Calculate 
Jj.(final  J)  and  from  Fig,  A-2  find  M^p.  Calculate  an  average  by 
means  of  the  following  formula. 

ml  *  °-03  (Mi.oZ  -  mlf!>  +  °-5  <mlo  +  mlf) 

where  M  is  the  Mach  number  at  the  end  of  the  grain,  and  the  subscripts 

"O"  and  *F"  again  refer  to  initial  and  final  conditions,  respectively. 

From  Fig.  A -3  find  the  value  of  P  /P  corresponding 

ti  c 

to  Mt  and  substitute  in  the  formula  for  C  , 

L  13 

o 

^  x  =  web 

*S  =  ScLx  ,  where  3  =  surface  area  (in2)  and  web  =  burning 

x  -  0 _  distance  (in).  For  cylindrical  grains, 

Grain  volume  (in3) 
in  web  (in) 


web 


number  as  a 
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c+  (ft/sec)  *•  characteristic  exhaust  velocity  =  ^r~— 

Do 

lUbf-scc/lbm)  -■  delivered  specific  impulse  =  J~- ^  (Ikf-flec) 
spa  '  '  11  charge  wt.  (lbm) 

,  h,  i  ,  ,  charge  wt.  (lbm)  X  1000 

m  (lbin/sec)  =  mass  discharge  rate  =  - ® - > - 5 — r-4 - 

'  t  (msec) 


C_  =  measured  thrust  coefficient  *=  I  flbf-sec/lbm)  X  C  (lbm/lbf-sec) 

spa  Dd- 

C_  =  corrected  thrust  coefficient  =  F?ooo  (lbf-sec/lbm)  X  C_  (lbm/lbf-sec) 
*  1000 

F?mo  =  specific  impulse  corrected  to  1000  psia  chamber  pressure,  optimum 
expansion  ratio  at  sea  level  atmospheric  pressure  (14.7psia),  and 
0*  nozzle  exit  divergence  angle. 


. ,  rvfi>‘ 

.p^Cp(ta.O 

Cp  (std)  is  a  function  only  of  specific  heat  ratio  {y)  and  is 
taken  from  a  table  of  such  values. 

Cp  (test)  is  calculated  from  the  relationship  > 


CF  (test)  =  XCF  (vac)  +  (1  -  X) 


XR 


where  C  (vac)  and  P  /P  arc  functions  only  of  y  and  XR 
JP  6  a 

(nozzle  expansion  ratio)  and  are  ob  ined  from  tables  of  such  values; 

P  ^  is  the  average  atmospheric  pressure  at  which  the  firing  is  made 
atm 

(14.4  psia  at  Huntsville,  Alabama);  P_^  is  the  actual  motor  average  pressure 
over  the  action  time,  and  X.  i  •  the  nozzle  divergence  angle  correction 

t 

factor  defined  by 


X.  =  Vz  0  +  cos  p) 


where  p  is  the  divergence  half  -angle  at  the  nozzle  exit. 
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=  Ffoo,,  when  the  following  additional  conditions  have  been  met: 

(a)  the  trace  is  neutral  within  the  limits 

0.90  £  £-  £  1.10 
b 

over  the  equilibrium  portion  of  the  trace, 

(b)  800  £  ?a  £  1200 

(c)  the  tail-off  portion  of  the  trace  is  limited  by  the  following: 

t,  S  0.87  t 
b  ap 


tfooo 


4.  Exceptions  for  Micro-Motors 

Some  definitions  must  be  revised  to  provide  additional 
flexibility  for  application  to  micro-motors.  The  trace  time  ^points 
are  still  applicable,  but,  the  location  of  these  points  will  depend  on  the 
shape  of  the  trace  to  be  reduced.  Several  possible  forms  of  pressure*-* 
time  traces  which  may  be  obtained  from  micro-motors  are  shown  in 
Fig.  A -4. 

The  exceptions  for  micro-motors  are  outlined  below; 
other  data-reduction  procedures  remain  unchanged. 


/Pbdt 


0.95 


«  I  #  «  Pjooe  when  the  following  additional  conditions  have  been  met: 

(a)  900  a  p  a  lioo 

a 

(b)  8  a  xr  a  1 1 

(c)  the  record  is  obtained  from  a  6CS-11.4  static-test  motor. 


I  1 

Bps 


lJANAF  Solid  Propellant  Rocket  Static  Test  Panel,  'Recommended 
Procedures  for  the  Measurement  of  Specific  Impulse  of  Solid  Propellants,  * 
SPLA  Publication  No.  SPSTP/lOB,  March  1962. 


A-l  1 


I 

I 

I 

I 

I 

I 

I 


Parameter 


Ignition  delay  time  (msec) 

(msec)  =  burning  time 

t^pCmsee)  =  action  time 
t  (msec)  *  correction  time 

P^psia)  -  burning-time-average  pressure 


P  (paia)  s  action-time-average  pressure 

fit 


3?  (psia)  *  equilibrium  pressure 
eq 


Integral  ratio 


Rise-time  factor 


Decay-time  factor 


Definition  and  Comments 
t2  -  tj 

ij  -  14  (use  for  burning  rate  and 

P,  calculations) 
o 

t8  -  tg  (use  for  P  calculations) 
& 

t8  -  tg  (use  for  beat  loss 
correction) 

+  14.4  (use  for  P-r 
17  ~  **  correlations) 

+  14.4  (use  for  C_ 

8  *  9  test  calculation) 

*j&~-  +  14.4  (use  for  P-K 
7  ”  ^  correlations) 

‘fofs'ft  (indication  of  "square* 

2  trace) 

(indication  of  poor 
17-12  ignition) 

t!  "~‘  I'7  (Indication  of  poor 
*  6  ignition  and/or  off- 
center  mandrel) 


i 


1 

i%3 

J 


•S 

% 


I 

A 

I 

1 

1 

I 

I 

i 

1 
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Appendix.  B 

Calculation  of  Velocity  Lag  In  Two -Phase  Flow 

A  digital-computer  program  written  in  FORTRAN  for 
calculating  particle-velocity  lags  in  two -phase  flow  through  nozzles  of 
specified  geometry  was  developed  using  a  simple  mathematical  model 
to  approximate  the  flow  process.  The  mathematical  model  1b  described 
and  procedures  for  using  the  computer  program  are  outlined.  A 
listing  of  the  program  and  an  example  problem  are  given  in  Appendix  C. 
Source  decko  for  the  program  are  available  on  request. 

1.  Mathematical  Model 

The  mathematical  model  describing  two-phase  flow  through 
nozzles  was  obtained  from  a  simplification  of  a  previous  study  at  this 
Division1.  In  this  simple  model,  the  combustion  and  expansion  processes 
were  assumed  to  be  adiabatic.  Dynamic  non -equilibrium  effects  were 
determined  by  treating  the  particle-velocity  lag  as  a  perturbation  of  a 
system  in  which  the  particles  were  assumed  to  be  gaB~lik@.  The  one- 
dimorisional,  icontropic  expansion  of  the  gae-like  system  was  first 
determined  through  the  use  of  an  effective  specific-heat  ratio  andan 
effective  gas  molecular  weight  based  on  a  chemical  composition  frozen 
at  the  combustion-chamber  conditions.  The  velocity  of  the  condonsed- 
phase  particles,  which  were  assumed  to  be  spherical,  was  then  approximated 
by  calculating  the  velocity  of  a  single  particle  introduced  into  the  gas- 
like  system.  The  coupled  effect  of  the  particle  lag  on  the  flow  of  the 
gas  -like  system  was  ignored. 

The  specific-hcat  ratio  of  the  gas -like  system  is  given2 


i 

'Rohm  &  Haas  Company,  Quarterly  Progress  Report  on  Engineering  Research, 
No.  P-60-17,  December  I960. 

zUrJted  Technology  Center,  ’’High  Performance  Nozzles  for  Solid  Propellant 
Rocket  Motors,1  TTR -7.026,  November  196't. 


i^jSkiAisU't 
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by 


V  (1  +  zr) 
i  +  v  zr  * 
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whore  Is  the  specific  heat  ratio  of  the  gas  phase  in  the  chamber, 
z  is  the  ratio  of  the  mass  flow  rate  of  the  particles  to  that  of  the  gas, 
and  r  is  the  ratio  of  the  particle  specific  heat  at  constant  pressure 
to  that  of  the  gas.  The  effective  molecular  weight  of  the  gajs-like  system 
is  given  by 


M  =  M  (1  +  z)  , 

0  g 

where  is  the  molecular  weight  of  the  gaB  in  the  chamber. 

Motion  of  the  particle  was  assumed  to  be  governed  by 
the  differential  equation 
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r  - 

du 

p 

- 

?  °ppP] 

U  = 

p  dx 

D 

Z 

Pg 

Mm* 

_ a 

D* 

P 

<*>  <u„  -  u  )|u  •  u  J  B-3 

g  P  B  P 

where  u  is  the  particle  velocity,  u  is  the  velocity  of  the  gas -like 
P  g 

system,  p  is  the  pressure  of  th  a  gas -like  system,  x  is  the  distance 
S 

along  the  nozzle  center  line,  is  the  particle  diameter,  is  tho 
particle  density,  is  the  density  of  the  gaB -like  system,  is  the  drag 
coefficient,  and  <{>  is  the  Cunningham  correction  for  molecular  mean-free  - 
path  effects.  The  positive  direction  is  assigned  to  the  direction  of  flow 
of  the  gas -like  system. 

A  least-squares  fit  to  the  drag  coefficient  tabulation  for 
spheres  given  in  Perry1  was  used  in  the  computer  program.  This  expression  is 


lJ.  H.  Perry,  ed. ,  "Chemical  Engineer's  Handbook,  *  3rd  ed.  ,  McGraw- 
Hill,  New  York  (1950). 
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inC  =  3.271  -  0.8893  (inRe)  +  0.03417  (inRe)2  +  0.001443  (lnRe)3f,  B-4 

0.1  £  Reg  2000 


where  the  particle  Reynolds  number  is  given  by 

D  (u  -  u  )  p 
Re  =  _El£ - £L_iL 

Hg 
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The  viscosity  of  the  gas -like  system,  p.^,  was  determined 
as  a  function  of  temperature  using  a  mole  fraction  weighted  average  of 
the  viscosities  of  the  major  gaseous  constituents  of  the  combustion  products 
of  a  typical  C-H-O-N-Cl-Al  propellant1.  The  viscosity  so  determined 
is  represented  by  the  polynomial 


in(p  X  105)  *  0.3757  +  0.8629in[r^]  -  0.05446^  , 
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1000*R  2  T  <  6500‘R, 

_  s 

where  viscosity  is  in  lbm-ft"l~sec  *  and  T  is  the  gas  temperature  in  *R. 


-l 


g 

The  Cunningham  correction  is  given  by*'* 


4>  '=  - r - r - 1 ,  B-7 

1  +  [  1.644'+  0.552  exp  (-0.656D  /\)j 

P  P 


where  the  units  for  the  particle  diameter,  D" .»  and  the  mean  free  path 
of  the  gas -like  system,  \(  must  be  identical. 

The  mean  free  path  of  the  gas -like  system  was  based 
on  a  mole  fraction  weighted  average  of  the  gaseous-phase  constituents. 


'Rohm  fct  Haas  Company,  Quarterly  Progress  Report  on  Engineering  Research, 
P-60-17,  December  I960. 

2J.  H.  Perry,  ed.  ,  “Chemical  Engineer  *s  Handbook,  *  3rd  ed, ,  McGraw- 
Hill,  New  York  (1950). 

3C,  E.  I  apple,  "Fluid  and  Particle  Mechanics,"  1st  ed. ,  University  of 
Delaware,  Newark,  Delaware  (1954). 
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The  mean  free  path  so  determined  ia  given  by1'2 

T 

X.  =  1.025  X  10"6  — ,  B-8 

'  8 

where  mean  free  path  ia  in  ft,  temperature  T  is  in  *R,  and  pressure 

8 

P  is  in  lbf-ff2. 

8 

The  nozzle  geometry  employed  in  this  analysis  is  shown 
in  Fig.  B-l.  The  nozzle  is  composed  of  a  conical  converging  section 
having  an  inlet  diameter  Dq  and  a.  half-angle  a,  a  rounded  throat  section 
having  a  radius  of  curvature  R  and  a  throat  diameter  Dfc,  and  a  conical 
diverging  section  having  a  half-angle  (3  and  an  exit  diameter  D^.  The 
coordinates  of  several  reference  points  in  the  nozzle  are  shown  in 
parentheses.  Similar  nomenclature  ia  used  to  describe  the  nozzle 
in  the  computer  program. 


Fig.  B-l  Nozzle  geometry. 


'Sir  James  Jeans,  'An  Introduction  to  the  Kinetic  Theory  of  Gases," 
University  Press,  Cambridge  (1948). 

2J.  O.  Hirschfelder,  C.  F,  Curtiss,  and  R.  B.  Bird,  "Molecular  Theory 
of  Gases  and  I.iquids,  "  Wiley,  New  York  (1954), 
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2.  Computer  Program 

Solutions  to  Eq.  B-3  were  obtained  by  numerical  integration 
using  a  Runge-Kutta  method.  The  apace  iir  the  conical  converging  section 
was  divided  into  unequal  increments  using  a  weighted  Mach  number  increment 
to  provide  smooth  starting  for  the  integration.  The  space  in  the  rotmded 
converging  section  was  similarly  divided.  The  space  in  the  rounded 
diverging  section  was  divided  into  unequal  increments  such  that  the  area 
ratio  at  each  successive  interval  increased  by  0.05.  The  space  in  the 
conical  diverging  section  was  similarly  divided,  except  that  the  area 
ratio  at  each  successive  interval  increased  by  0.25. 

The  units  of  the  input  data  for  the  computer  program  are 
those  commonly  employed  in  internal  ballistics;  the  computer  program 
converts  the  data  to  an  internally  consistent  set  for  calculation.  The 
data  required  by  the  computer  program  are  described  below* 

_ Mnemonic  Symbol  _ Description _ 


DZERO 

Nozzle  inlet  diameter,  in. 

R 

Throat  radius  of  curvature,  in. 

ALPHA 

Converging  section  half-angle,  degrees 

BETA 

Diverging  section  half-angle,  degrees 

CAMMA 

Effective  specific  heat  ratio  (Eq.  B-l) 

TZERO 

Chamber  temperature,  *K 

PZERO 

Chamber  pressure,  psia 

XMW 

Effective  molecular  weight  (Eq.  B-2) 

CJ.NTV 

Number  of  integration  intervals  in 
conical  converging  section  (usually  20) 

XR 

Expansion  ratio  at  exit  plane 

DP 

Particle  diameter,  microns 

DT 

Throat  diameter,  in. 
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The  output  format  consists  of  a  heading,  a  listing  of 
the  data,  and  a  tabulation  giving  the  particle  velocity  in  ft-scc"1,  the 
local  expansion  ratio,  and  the  relative  velocity  lag  (u  -  Up)/Ug  as  a 
function  of  the  distance  in  inches  from  the  nozzle  inlet. 


Appendix  C 

FORTRAN  Program  for  the  Calculation  of  Velocity  Lag 


00000  * _ VELOCITY  LAG  IN  TWO  PHASE  PLOW _ 

ooooo  c 

-0-0  0-00 C _ P  BJS.QftAJSJUEJLBl—L* EL,_SX0.NEC.tP.H£&. _ _  _ 

00000  C  R0HH  and  HAAS  COMPANY 

_0_0  00  0  c  HUNTSVILLE.  ALABAMA _ 

ooooo  c 

-QJPOOO _ C _ _ _ _ _ __ _ 

00001  ~  COMMON  ET  A  l »  ETA2*  ETA3  *  £TA4,““ETA5.  ETA6.  ETA7. 

00001  2  2TA6,  T7ERO »_RHO_Z£  P_N£_N1.  ROOT ,  j  TESX-f— UQAS  . 

00001  2  XMACH,  HMACH,  YMACH.  xi,  X2 «  Y 1 ,  Y2.  XXR.  YXR.XUP 

00001  1 _ WE A D?>DZERQ,R, ALPHA* BETA, GAMMA . T ZE RQ , PZER 0  > _ 

00001  2XMn.C!NTV,XR.DP.DT 

00034  2  _ P 0 RMAT  <5E1*. 8 ) _  _  _  _ 

00043  Sl'a  DZERO 

00047 _ S2_s_R _ _  _ _ 

00053  “S3  *  ALPHA 

00057  S4  =  BETA 

00063  s£  =  GAMMA  ' 

00067 _ S6  *_TZER0 _ ' _ 

0  00  73  S7"~*  PZERO"  “  . '  " 

00077 _  S8  *  XMW 

00103  DZERO  a  D ZERO/12 .0  " 

00111  R  =  R/12.0 

Torn  iLPM A  “  alPhA /57'. 295 76 

00125  BETA  *  8ETA/57. 29578 

00133“  '  "  "TZE'ft0‘  “  fZE'R0*r.  ii 


00141  PZERO  =  PZERO«144.0 

D  rri n  z~  p  pint*  ,~s f.s  2  rs3;sx;s5  rs67 sr/s's 

00172  4  PORMAK 1M120X25H  VELOCITY  LAG  CALCULATION  /// 

1T0T72  2  ~^-ro  2?R<r  =  Fxr.'?  7”  6)f5n“R  =  P  9 , 5  ,  2X94  aIpVTa  =  f9.5“, 

00172  2  4X8H  BETA  *  P  9 . 5  /  9H  GAMMA  =  P9.5r2X9H  TZERO  5  T9. 2. 

0 0 172 - 2“7xWnPZER0 "  = ' 79  ;2T«tf SH'Fir'TTOTT^ - 

00356  S9  »  DT 

00362""' .  S10  V'DP  -  - -  - 

Q0366  DT  «  PT/12.0 

wtta - \)p  »  pyy3T7FT7TDwn (nrm - 

00*04  PRINT7.CINTV,XR,S10.S9 

00417“  ’7  . rURFTA7{T)f8HCIWTV'  FTOri"'-  6X5MXR'a  T9.5  T^VSffDP  ’  =  f9,5,“ 


00417  27X5HOT  t  r9.5  //) 

00470  C  PEPrORM' PPECOMPUTATI ONS 
00470  RHOP  *  3  .  Q7*62 , 4 


TJTT4T6 - 

0050? 

00520 

00534 

0055? 

00572 

"T  x  37TTT5V77 - - 

RHOZf  *  <PZERO*XKW)/(TZERO*1546.0) 

ETAl  =  2.0/CGAhMA  *1.0) 

ETA?  *  (GAMMA  -  1.0)/CGAMMA  ♦  1.0) 

ETA'3  *TGAMHA'VT.O)/(2.0  #(GAMMA  -  l.0)> 

ETA4  *  ETA3  -  1.0 

-- 

*00  503" 

mv  *  fG7(fl>n"  -  T7tn7275 

00610 

fc  T  A  6  «  1 . 0 / ( GAMMA  -  1.0) 

00624 

E  T  A  7  =  (1546.0*32. 1 7/ XM w ) *2 . 0 *G A MMA • E T A  6 

00642 

ETA*  »  3.0/ (4  .  n*OP*9HOP > 

00660’  e 

■  calculate  nozzle  parameters . *~  ’"  . . . 

00660 

CALPA  s  COS(ALPHA) 

00666 

sum  =  sjNULt’M.i 

00674 

TALPA  S  SALPA/CALPA 

4 


c-z 


TO 7C?  CBET  a  =  COS ( bS  T  A  J 

00710  _  SBETA  =  SIN(BETA) _ 

‘00716  .  TrETa  =  SBET* /CBET A  . . 

007?4  A  s  ( DZEPO  -  DT  ♦  2.0  *  R  *  <1,0_/  CALF*  -  1.0))/ 

00724  2  12.0  *  TALFAJ  . 


00764 

inmr- 

01005 

oioi6 ; 
01026 
01037  ” 
01055 
TJTTJ67 
01067 
01113 
01131 
01155 
01173 

trrm— 

01217 

01225 

01233 

01237 


fi  J  R  ♦  DT  /  2.0 

E  A  -  R  •  SAL?  A 
F  *  B_  -  R  *  CALF* 

V.:  G'  *  A*P  *  SSETA 

M  *  H  C8ETA 

qT/2 , 0  V*SORT  OlR ) 

XK  =  G  ♦  <  XL  -  H)/TB€Ta 

T - CTArcrOTA  T  E‘  “  a  R  ET'R  i  T 1 0  S’  *  — 

AREAO  s  IP1/4 .0>«D2ER0*»2 

. ARF'AF  =  pj  #F*»2 

AREA  T  =  ( P 1 /* . 0 ) *DT • • 2 
ARE  Ah  =  Pl*rS«»2 
A«EALsPl*XL**2 

- YR(T"=  ARE AU7VREAT - 

XRF  =  ARFAF/AREAT 

XRh  =  AREAH/AREAT -  ' 

XRL  -  >» 

GO  TO  45 


01241  C  SET  UP  PRINT  FOR  ERROR  WHEN  MACH  SUBROUTINE  ITERATION 

01241  counter  is* exceeded.*' - - 

01241  ft  PR | NT9 

01244  9  FORMaT  <  2Bm  ITERATION  COUNTER’  E’XcETDt’trT . . .  — ' 

01267  GO  TO  44 

01271  C  .  BEGIN  CALCULATIONS  IN  CONVERGING  "SECTION  W  N022Le - 

01271  45  ITEST  =  1 

GT775  . TJUESS'  i  0.01 -  - -  ‘ 

01301  CALL  naCH(GUESS.  XP0) 

01306  *  rr{ROOT-l. )  4ft .  46.46 - - - - - 

01316  46  PR1NT47 

01321  *7  FORMAT  ( TX17HC0NVERGENCF  ERROR  -  '  '  . 

01340  GO  TO  44 

Trt3*?— *8 - IF  CTTESTTfl  ,  8  ,  V3 - - - - 

01350  10  OMACH  b  ROOT 

0135* —  call  MACHfGUESs,  xnn .  . . .  . .  —  -  ’  - 

01361  IF (R0OT-1. 0)49,46.46 

01371  49  IF ( 1  TEST ) ft . 8  » 11 

013/7  11  FMiCH  =  ROOT 

TT1'4C3— t  CTCCUL*TF**»TE  I  GHTED  ’I  VCPEtfETIT-?!  2^  ITT  "CCNvFRG TNG"  SECT  TON’  " 

0  1 4 1]  3  0ELXbCINTV«(CINTV»1.)/2. 

01415  .  •  0ELM1*<FMACH-0MACh>/DELX  . . 

01425  D6LM2«<1.-FHACH)/DELX 

01435  C  PRFFrX  X  ‘  DENOTES  PRESENT  PARTICLE  LOCATION,  PRFFUY"* 
01435  C  OENOTES  INCREMENTED  LOCATION.  PREFIX  H  DENOTES  MID 

~0T4  35  C”  POTNT  'Of  X-Y'  INTERVAt ""SUrmf  T"ANTT  2~'C0RRFSPOTr_TO - 

01435  C  X  ANO  Y  POSITIONS,  RESPECTIVELY. 

01435  (T - SE1  DATA  FOR 'INITIAL  LOCATION  ’ 

01435  XI  s  0.0 

urrAl - YT  -  DZER0/2'.Q -  - - - - - 

01447  XXR  :  XRO 

“ott53 - mrrr  •  “<jMrcw - - 

01457  MM  S  1 


% 


% 

< s> 


°n 
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_ TE»M  =  1,0  »  ETA5  •  OHACH  ••  2 

Tlt03  TEnR  !  T2eRo/TeR H 

01511 _ UG  a  S  «  SORT(ETa7.(TZER0-TEHP>) 

01527  DELTA  i  0.001  "' 

01533 _ XJ°  i  0 . 999*UG A  S _ 

01551  PR1NT12. XI. XUP.XXR, DELTA 

01555  12 _  _FORHAT(2OxlHx,lOx2HuP,10x2HxR,9x5HOELTA/18x5H(]N. 

01555  2.3.F12.1,  F12.3.  F12.5  > 

0163J) _ C _ BEGIN  DO  LOOP  FOR  CONVERGING  SECTION 

01630  CVcTnTV - - - 

01635  13  DO  25  I«1,L 

OliAl  X!  =  ]  - - 

01656 _ GO  TO  (15,15)  .  HN 

^16?5  14  YMACH  «  XHACH  ♦  DELH1»XI 

0166_4_  _  GO  TO  16 

SlSSi  15"!  fHACH  »'  XMACfH  ♦  DELH54XI 

01676  16  AREA2  «  ( AREAT/YH ACH) • ( ET A1*ETA2»YHACH» *2) 

tots - rrrmr- - - 

01735  Y2  *  SORT ( AREA2/PI ) 

0175  i  G o~To  TIT, IS)  ,HH 

01757  17  »;■(!■  ( Y 1  -  Y2)/TALFA _ 


(  FPS ) / 17XF6 


02061 
0207"S 
02101  20 
imT* 
02116  21 


X2iA-S0RTlA8S(R»»2-(6-Y2)««2) ) 
YxR  «  A9EA2/AREAT  ' 

HX  «  ( X2  -  XD/2,0  a  XI _ 

“00"  to  (20.21)  .MM 

HY  «  Y1  -  T ALFA* ( HX  -  XI) 

GO"  75  5?  - - 

HY  *  B  -  S0RT<R«»2-(A-HX)«52) 
"RXR  ■  'dn""HV»»2 WlftEAT 
GUESS*XHACH 

■"nirriim  - 

IF(ROOT-l. 0)50.46, 46 


02237  23  HMACH  *  ROOT 

07257-0  USE  SUBW8UTTRF  "TO"  CS'LLTOTE  UP  AT 

02243  CALL  FINAL 

0225  6  75  CONTTITOF - - — - — — 

02253  C  END  OF  CONVERGING  SECTION  DO  LOOP 

S225J— - ir(HK--7T25725r76 - 

02263  25  MM  ■  HH  ♦  1 


02273  C  TWO  AREA  INCREMENT  SIZES  ARE  USEO  IN  DIVERGING 

0227T  c - secTioN-or-nmcr"" — - — 

0227326  XN  •  ( XRH  -  1.01/0.25  ♦  1.0 


02332  HH  ■  1 

"02336  2 T - HD"  37~T*  l.ffN - 

02342  00  TO(26. 29.30). HH 


02367  30  Y2  «  SORT ( VXR«ARE AT/P  I ) 


1 


C-4 


Q2A  '  6 
02413 


CU4SS~"*H»Ch  »  SiST" 
CALL  HACHIGUESS.  YXR> 


H717JS - TFT1.  0-ROOTT5T,  46,  «- 

02430  51  I F ( I TESTI8.8.31 _ 

T2TJ* — SI - W'  9OTT 

024  42  1 F  t  Y  a  -  H132.32.33 

0-24 52 — S2 - r5‘*  1  '•  ‘SORI  (h*»2  ~ 

02523  00  TO  34 

TTZ525 — S3 - 22'* 


02551 

02561 

35 

— TFTHr  ^'  G  )  35 ,  « .  ab 

MY  >  B  -  SORTtR  ••2-(*-HX)»*2) 

0  V  6  i  -i — 

"  BU  10  il  ,  „ 

Q2635 

36 

HY  «  H  ♦  THETA  •  <HX  -  51  _ 

“0Z6<T~ 

02667 

T7 — 

GUFSS*(XMACH*YMACHW2.0  _ _ 

TT2577 

- - 

— CTtt'WAUHlUUFSS'T-  HTTfT' 

02704 

IFll.0-ROOT152.46.46 

JJTTPT 

”57 

- 1*  l I  I hb 1  )0«y *30 

02722 

30 

MRACH  »  ROOT 

”07776  final 

02736 

»f  <Mm-2MOa60,^  .  _ 

U  2  7  4  6" 

— <U 

- ffl-T  FIR  *  1 

02794 

NN  «  XNN 

027  sr 

- - 

— go  io ~'n 

02763 

60 

MM ■ HH  * 1 

0277T" 

— YTX*r* - - 

02775 

NN  ■  1 

osinrr 

- BTTTCT77 

03003 

44 

00  To  1  - 

'  (T3TT0V 

* 

— nc  _ 

Final  Subroutine 


00000 

• 

FJNAL  SUBROUTINE 

00001 

00001 

— 

Subroutine  final 

COMMON  ETAl.  6TA2,  £TA3,  ET  A  4 .  FTA5*  ETA6,  ETA7.  eTAB. 

0000* 

00001 

JTfSRO,  RHOZE.  DP,  FNCNi.  ROOT.  1  T  8  S  T  »  UGAS.  XMACH, 

2  HHACH,  YMACH,  Xl»  X2,  Yl,  Y2,  XXR,  YXR,  XUP 

00001 

c 

Begin  rUnge  -  kutta  Integra t ion  procedure 

00001 

Call  fncn(xmach.  xup> 

00006 

Cl  *  ( *  2  -  XI ) «  F  NCN 1 

00016 

ARG1*C1*XUP**2 

00034 

IF (ARG1)200»203.203 

0004? 

203 

HUP  a  SORT (C1*XUP**2 ) 

00066 

CALL  FNCNIHMACH.  HUP) 

00073 

C2  s  ( X2  -  XI' *FNCN1 

00103 

ARG2=4.0«C2-2.0*C1*XUP«*2 

00137 

IF( ARG2J200 .204 . 204 

~t)  0 145 

9TTT 

yup1=sort  <27o*(2 .  o»cTFci  )«xup»*'2'> 

00203 

call  fncncyhach,  YUP1) 

00210 

C3  *  1X2  -  X 1  )  *F  NCN 1 

00220 

»RG3*XUP*«2*(C1*4.0*C2*C3 1/3.0 

0025? 

IF (ARG3 1 200 . 205 , 205 

00260 

20? 

YUP*  SORT ( XUP* *2* (Cl*4 . 0*C2>C3>/3.0) 

00320 

C 

END  OF  INTEGRATION  PROCEDURE.  check  for  up  GREATER  than  uc 

00320 

I F  <  UG AS  -  YUPJ200, 200 , 201 

00330 

200 

YUP  =  0.99V*  UGAS 

00336 

201 

DELTA  «  (UGAS  -  YUP )  /UGAS 

TSJ45- 

xp  *  x2«T576 

00354 

C 

PRINT  FROM  Y  LINE.  NOTE  X  UNIT  CHANGE  FOR  PRINT 

00354 

PR  INT202.XP.  YIJP, YXR,  DELTA 

00367 

202 

F OHM aT ( 17XF6 , 3 .  F12.1,  F12.3,  F12.5] 

obno 

C 

bATA  Transfer  fOoh  v  line  To  x  line 

00  4  3  0 

XI  <  X2 

06414 

Tl  «  Y2  ...  -  .  . . 

00420 

XUP  s  YUP 

■604?* 

XxR  *  yxT 

00430 

XMACH  *  YMACH 

TO  434 

return 

00440 

• 

END 

Best*’'*16'8 


Cop' 


Function  Subroutine 


00000 

• 

FUNCTION  SUBROUTINE 

“Bwirr" 

00001 

Subroutine  fncnixmock.  up) 

COMMON  ETll,  ETA2.  ET  A  3  »  ETA4,  ETAS.  El A  6  »  E  TA7  . 

00001 

00001 

^  rrifl,  t  zero .  rhoztt^ d  p,fn c  ni ,  root. 

2  XMACH.  MMACH.  YMACH.  XI#  x2.  Yl.  Y2. 

1  TEST  ,  UGAS. 

XXR »  YXR .  XUP 

00001 

00001 

C 

EVALUATE  RUNGE  -  KUTTA  FUNCTION  IN  EQ 
TERM  *  1.0  ♦  ET A 5  •  X MOCK  *»2 

3 

00030 

TEhT1  ■  TZEHO/TERM 

00036 

RHO  *  RH0ZE/TERM#*ETA6 

00063 

00071 

UGAS  *  SQRT(ETA7»(TZER0  -  TEMP)) 
T6RHl*Et-OG<TEMP/100  0  .  ) 

“irofoB- 

0Q105 

XMU  0.8629  •  tERMl  -  'f. 

2  TERM1**2>/100000 .0 

65446. 

0QH6 

00162 

XLAHB  =  2.052/<RmO*1.0E08) 

RATIO  *  XLAMB/DP 

oTim 

P"H i  t~X .  (T/ a  .T  ♦  Ra TTS •  i  l .  6~4  4  ♦  0 . 5S2 

.  EXP(-0.456 

00170 

2  /RATIO)  )  ) 

0  6223 
00250 

RE  *  ABS< DP*( UGAS  -  UP ) •RHO / X MU)  ~ 

RELOG*ELOG < RE  ) 

T5 T-T* 

T)R»G  *  EXP<X271  -  07eT?3  •  RELOG  •  0 

■fj417  • 

00256 

2  RELOG  **2  ♦  0.001443  *  RElOG  •**) 

0"QT3T“ 

FUG  Ml  *  TJBac  •  RtJfl  •TT |5*S  -  UR! 

•  PHI  ♦ 

0  0334 

2  ABS ( UGAS  -  UP) 

00402  “ 
00406 

• 

RETURN 

END 
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00000 

OOOrti 

OOOCl 

00001 

00001 

TJUMT" 

OUOOl 
0  U  0  0  *? 
00034 
00043 
00043 

ovm 

00133 

00156 

00162 

coin 

00177 

-mrr 

00207 

0021' 

002lo 


Mach  Subroutine 


*  _ MACH  SUBROUTINE 

"SUBRuUTINE  ‘H*CH<«OOTi.  RATIO) 

COMMON  E  T* 1 ,  ETA2,  ET  A3 »  FT  A  4 ,  ETA5.  6TA6,  ETA7, 

2  ETAS »  T2ERO.  RHOZE.DP#  fNCNl.  HOOT,  1  TEST ,  UGAS, 

2  XMACH  #  HMACH.  YMACH.  XI,  X?,  Yl ,  Y2,  XXR,  YXH ,  XUP 

*r  ”  "  CALCaiirF'HXCR  NOKBEH  . . . . 

K  *  1 

l Uu  Term  1  (f-TAl  ♦  6TA2  *  R00T1  •  •?) 

TERH1  s  TERH*??Ta3 
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Appendix  O 

Velocity  Lag  Curvea  for  RH-P-1  12  Propellant 


Curve#  of  A  versus  nozzle  throat  diameter  for  particle 
diameter  (N)  values  of  0.5,  1.0,  1.5,  and  2. Op  and  nozzle  expansion 
ratio  (XR)  values  of  6,  8,  10,  and  12  are  shown  in  Figs,  D-l  through 
D-6  for  the  Rohm  h  Haas  .730.50,  2C1.S,  and  6C5  motor  nosales; 
for  the  USAF  BATES  motor  nozzles;  and  for  the  Aerojet  10KS-2S00 
motor  nozzles.  These  curves  were  generated  using  iho  FORTRAN?  program 
described  in  Appendices  B  and  C.  The  thermodynamic  input  parameters 
used  were  those  determined  for  RH-P-1 12  propellant  (QAMMA  ■  1,141, 
TZERO  *  3389‘K,  XMW  «  27.452).  The  parameter  PZERO  ■  1000  psla 
was  used  in  all  cases.  The  nozzle  configuration  parametere  used  are 
listed  in  Table  D-I. 


Table  D-l 

Nozzle  Configuration  Parameters 


Motor 

.75C.50 

2C1.5 

6C3 

BATES 

1  OKS -25  00* 

DZERO 

0.75 

2.00 

6.00 

11.73 

7.125 

R 

0.4375 

1.167 

3.50 

3.00 

4.0 

ALPHA 

45 

45 

45 

45 

30 

BETA 

15 

15 

15 

15 

15 

‘Estimated  from  non-dimensioned  drawings. 


75C.50  motor 


NOZZLE  THROAT  DIAMETER 


Fig  D-2  Velocity  log  curves  for  .75C.50  motor  nozzles. 
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NOZZLE  THROAT 


08  .10  .12  .14  .16 

FRACTIONAL  VELOCITY  LAG 

Fig.  D-3  Velocity  lag  curves  for  20.5  motor  nozzles. 


i 


NOZZLE  THROAT  DIAMETER  -  in. 


4.0 

3.5 

3.0 

£5 

2.0 

1.5 

1.0 


.01  .02  .03  .04  .05  06  .07  .08  .0 


FRACTIONAL  VELOCITY  LAG  -  Z 


Fig.  D-4  Velocity  lag  curves  for  6C5  motor  nozzles. 


NOZZLE  THROAT  DIAMETER 


Fig.  D-5  Velocity  lag  curves  for  BATES  motor  nozzles 
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RACTIONAL  VELOCITY  LAG  -  A 


NOZZLE  THROAT  DIAMETER 


Fig  D-6  Velocity  lag  curves  for  I0KS-2500  motor  nozzles. 
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Example  of  Scaling  Factor  Determination  and  Specific  Impulse  Prediction 


1 .  Scaling  Factor  Determination 

Assume  that  the  data  shown  in  Table  E-l  haB  been 
obtained  for  a  certain  propellant  through  static  teats  of  micro-motors 
at  =  1  000  psia. 

Table  £-1 

Micro-Motor  Static  Teat  Data 


Motor 

Grain 

Lcmgth 

in 

K 

m 

in 

ri> 

iit/oct: 

XK 

t 

X 

aec 

A,2 

in2 

in 

lbm 

p0 

-  *  IQQB. 

246.4 

.75C.50- 1.5 

1.41 

173 

0.146 

0.60 

B.'J 

0.?,34 

1.045 

0.02  189 

.75C.50-3.5 

3.41 

171 

0.225 

0.60 

8.2 

0*237 

1.053 

0*05241 

252.2 

2C1.5-1 

3,88 

170 

0,386 

0.60 

a.o 

0.470 

10.42 

0.3306 

252,6 

"Obtained  i\»  outlined  in  Appendix  F, 


Assume  further  that  the  values  listed  in  Table  E-II 
have  been  obtained  from  thermochemical  calculations  for  the  propellant. 


Table  E  -II 

Thermochemical  Data 


Equilibrium 

Chamber 

Chamber 

Chamber 

T  ,  *K  T  ,  »K 

c  e 

mols  gas/100  g 

mols  AI2O3/I 

3629  2224 

1 .20 

3.65818 

0.21239 

The  following  parameters  necessary  for  the  calculation 
of  fractional  velocity  lag  (A)  and  for  scaling-factor  determination  may 


E -2 


now  be  determined: 


*  = 1  -  t£  * 1  -ffff  c  °-38716 

c 


(as  sinning  C  =  C  ) 


p  p 

1  c  1  e 


iVi 

x  b  -i - —  =  1*0196  (mole  AlaOs/lOO  g) 

m  +  m 
s  g 


x  0.21655  _  „ 

z  ■=  - -  s  - -  .  j1  =  0. 27640 

1  -  x  1  -  0.21653 


1.0196  (0.21239)  *  0.21655 


V1 


1,986  X  (mole  gaa/lOO 
100 


1.2  1(1.986  X  3.65818 

1.2  -  1  100 


■  0.4359 


PA12Os 


0.4359 

0.344 


J..2672 


GAMMA 


\  (1  +  Zy)  1.2  {l  +  0.2764  X  1.26721 


l  +  v  »r 
g 


1  +  1.2  X  0.2764  X  1.2672 


=  1.141 


XMW  =  - 1  I0°  /.on"  =  a  =  ji7‘336 

mole  gaa/lOO  g  3.65818 


The  preceeding  information,  together  with  the  nozzle- 
configuration  parameters  listed  in  Table  D-I  would  provide  the  required 
input  data  necessary  to  calculate  A  values  for  this  propellant  using  the 
FORTRAN  program  described  in  Appendices  B  and  C.  However,  the 
curves  of  Appendix  D  will  be  used  to  obtain  A  values  for  this  example. 

From  curves  D-l  through  D-3  Table  E-HI  may  be 


formed. 


Tabic  E-III 


E-3 


A  Value  from  Curves  of  Appendix  D 
.  A  Values 


Motor 

t 

XR 

N  0,5 

N  -  1.0 

N  r.  1.5 

N  s  2.0 

.75C. 50-1.5 

0.146 

8.5 

0.0722 

0. 1 58/t 

0.2356 

0.3025 

.75C.50-3.5 

0.225 

8.2 

0.0490 

0.1143 

0.1811 

0.2398 

2C1.5-4 

0.386 

8.0 

0,0298 

0,0738 

0.1228 

0.1703 

We  are  >ow  in  a  position  to  calculate  q  valuea  from 
Eq«  29  in  the  text,  thereby  generating  the  q-veraua-N  curves  necessary 
for  Mealing -factor  determination. 


1  -  xAjI 
-  1  -  xAjj 

i 

(F?ooo)t 

2.T  n 

1 >  -  *A,f  | 

»/|  (“/l ' 

™T0I 

\i  -  xAt|  ' 

* 

i  m  j‘\  m  1 

1  z  - 

Taking  motor  i  »>*  the  .75C.50-1.5  motor  and  motor  2 
as  the  .75C.50-3.5  motor, 


<) 


_ (252.2)*  . 

2  X  778,16  X  0,38716 
114  X  32.174 


1  -  0.21655  X  0.0490l 

1 1  -  0,2)655  X  0.0722 / 

(246.4)* 

I  -  0,21655  X  0.0490 

1  0.234  X  1.6451 

1>.237  X  1.653  " 
0.05241 

1  -  0.21655  X  0.0722 

\  0.02189  J 

63,604,84  -  (1.0102I)(60,712.96) _ 

0713005  I  U. 01 02 1)(  17.58474)  -  7.47492] 


1697.9 


For  N 


l.OfJ-, 


63,604.84  -  i 

f  1  -  0.21655  X  0.1143 

2 

(60,712.96) 

i  1  -  0.21655  X  0.1582| 

0.13005 

1  1 

-  0.21655  X  0.11431 

-j 

(17.58474)  -  (7.47492) 

I  - 

-  0.21655  X  0.1582 

q  a  1243.3 


For  N  =  1.5^, 


_  63.604.84  -  1,02501  X  60,712,96 _ 

q  “  0.13003  (1.02501  X  17.58474  -  7.47492) 

q  a  1001.1 


For  N  =  2.0|jl, 

f  -  63,604,84  -  1.02927  X  60,712.96 _ 

1  "  0.13005  (1.0292?  X  i?. 58474  -  7.47492) 


q  a  806.8 


We  now  take  motor  1  a8  the  .75C.50-1.5  motor  and 
motor  2  aa  the  2C1.5-4  motor. 


For  N  =  0.5^., 


(252.6)*  - 

1  -  0.21655  X  0.0298 
1  -  0.21655  X  0.0722 

2 

■  (60,712.96) 

0.13005 

1  -  0.21655  X  0.0298  ‘ 

j - 

(17.58474)  -1 

0.470  X  10.42  ' 

1  -  0.21 655  X  0.0722 

0.3306 

J 

-  63,806.76  -  1.01873.  X  60,712.96  _ 

q  "  0.13005  (1.01873  X  17.58474  -  14.8136?)  ~ 


For  N 


l.Ofi. 


6 3.806.76  -  1.03820  X  60,712.96 _ 

q  =  0.13005  (1.0382  X  17.58474  -  14.81367) 

q  -  172  9.9 
For  N  =  1.5(i, 

63.806.76  -  1.05214  X  60,712.96 _ 

q  "  0.13005  (1.05214  X  17.58474  -  14.81367) 

q  =  -149.6 


For  N  *»  2.0(i, 


63,806.76  -  1.06220  X  60,712.96 
q  =  0.13005  (1.06220  X  17.58474  -  14.81367) 

q  -  - 1  3  5  8.  U 

The  q  values  calculated  above  are  plotted  against 
the  respective  assumed  N  values  in  Fig.  E-l,  Two  curves  result, 
one  for  the  .75C.50-1.5  and  the  .750,50-3.5  motors,  and  one  for 
the  .75G.50-1.5  and  2C1.5-4  motors.  The  curves  intersect  at 
q  -  1180  Btu/sec-ft3  and  N  -  1.12p.  These  values  represent  the 
scaling  factors  for  this  particular  propellant. 

2.  Specific-Impulse  Prediction 

Assume  that  specific-impulse  predictions  are  desired 

for  the  propellant  in  the  6C5-11.4,  BATES,  and  1QKS-2500  motors. 

Estimates  made  for  the  parameters  D^,  XR,  t2 ,  A,  and  m  for  these 

motors  (P  »  1000  psia,  K  —  170,  F  -  0.60  in/aoc,  density  -  0.062 
b  m  b 

lbm/i  n3)  are  summarized  in  Table  E-IV. 


EFFECTIVE  PARTICLE  DIAMETER  (n)  -  fl 


curves 


for  scaling  foe  tor  determination 
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Table  E -IV 

Estimated  Balliatic  Parameters 


Motor 

S  ,  in* 
m 

3t,  in 

XR 

web,  in 

t  ,  sec 
z 

A,  in2 

Volume,  in3 

m,  lbm 

6C5-11.4. 

194.78 

1.21 

8.0 

0.50 

0.94 

85.19a 

97.19 

6.026 

BATES 

620.48 

2.16 

8.0 

1.875 

3.53 

345. 89a 

1163.4 

72.13 

10KS-2500 

9 1 5b 

2.62 

8.0 

1.7b 

3.20 

102. lb 

1570 

97.34 

^Obtained  as  outlined  in  Appendix  F. 
“Estimated  from  non-dimensioned  drawings. 


The  next  step  is  the  determination  of  A  values  for  the 
.75C.50-1.5,  6C5-11.4,  BATES,  and  10KS-2500  motors  at  the  previously- 
determined  effective  particle  diameter  of  N  =  1.12|a.  This  may  be 
accomplished  by  plotting  A  values  obtained  from  the  curves  of  Appendix  D 
(for  the  specific  nozzle  configurations)  versus  N,  as  shown  in  Fig.  E-2.. 
From  Fig.  E-2  the  following  A  values  are  read  at  N  =  4,12p: 


For  the  .75C. 50-1.5  motor. 

A 

=  0.1785 

For  the  6C5-11.4  motor. 

A 

=  0.030 

For  the  BATES  motor. 

A 

=  0.017 

For  the  10KS-2500  motor, 

A 

=  0.014. 

Taking  motor  1  as  the  .75C.50-1.5  motor  and  motor  2 
as  the  larger  motor,  we  may  now  calculate  predicted  specific  impulse  of 
the  propellant  in  the  larger  motor  from  Eq.  28,  for  N  =  1.12p  and 
q  =  1180  Btu/sec-ft2. 


:lnh\8  COWY 

3est  Avauao"3 


t  A]  1  ,  ,  t  > 

_z  (  ZJtj  q  z 


(po  \  Ah[  -  xA;  (1<  c  \2  +  — —  -  HUiL 

^  1000)2  -Vtrr^J  ^l5«o)i  +  14.4ge  1  j  J  H4gG 


For  the  6G5-11.4  Motor: 


.-V  rF^irsIrFlj^')1  [l246--**1  +  <O.I300SK1  moniv.sa-iv-i)  (o.i.ioo5)(i i«o) 


(Fiooo)z  =  V  1.06800  (60,712.96  +  153.459  X  17.58474)-  153.459  X  13.28884 


U'tooo/z  s  Vl, 06800  (63,411.500)  -  2039.292 


(F?ooo)2  »  256.3 


For  the  BATES  Motor: 


-  0.21655  X  0.017  \z  ,,  ,  J3.53  X  345.89 

L-w  6X1655  XOTim  (63,411.500).  153.459  X 


(F?ooo)e  -  y  1.07406  (63,411.500)  -  153.459  (16.92765) 


(Fiooo)i  -  256.0 


For  the  1  QKS-2500  Motor: 


163  4 ]  )  500)  -  163  459  I— — 

-  0.21655  X  O.l  7851  97.34 


E-10 


(F^),  B  ▼  1.07547  (63,411.500)  -  153.459  (3.35648) 

(^1000)2  =  Z60.2 

To  summarize,  based  on  previously  determined  scaling 
factors  of  q  =  1 180  Btu/Bec-ft2  and  N  =  1.12p,  specific-impulse  valueB 
of  256.3  would  be  predicted  for  the  6C5-11.4  motor,  256,0  for  the  BATES 
motor,  and  260.2  for  the  10KS-2500  motor. 


i 
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Appendix  F 

Heat-Loss-Sur face-Area  Determination 

As  an  aid  in  determining  the  heat-loss-surfacc-arca 
term  (A)  used  in  Eqs.  28  and  29,  Figs..  F-l  through  F-4  have  been 
generated.  These  curves  show  nozzle  surface  area  versus  nozzle 
throat  diameter  for  the  .75C.50,  2C1.5,  6C5,  and  BATES  motor  nozzles. 
Nozzle  surface  area  as  used  here  is  defined  as  the  number  of  square 
inches  exposed  to  hot  gases  in  the  nozzle  convergent  section  only, 
beginning  at  the  nozzle  entrance  and  ending  at  the  nozzle  throat.  Nozzle 
surface  area  downstream  of  the  nozzle  throat  is  ignored.  The  total 
surface  area  (A)  exposed  to  hot  gases  Is  the  sum  of  the  area  of  the  firing 
head,,  area  of  initial  motor-case  ends  (not  covered  by  propellant),  average 
motor -case  area  exposed  as  the  grain  ends  burn  back,  and  the  nozzle 
surface  area. 

As  an  example  of  the  use  of  the  curves,  consider  a 
6C5  motor  having  a  nozzle  throat  diameter  (D^)  of  1.210  inch,  a  motor-case 
length  (L.^)  of  11.40  inch,  and  an  initial  grain  length  (L^)  of  11.25  inch. 
From  Fig,  F-3,  nozzle  surface  area  =  44.66  in2,  firing  head  area  = 

28.274  in2,  area  of  motor  case  ends  =  18.850  (11.40  -  11.25)  =  2.828  in2, 
and  average  area  exposed  by  burning  grain  ends  =  9.425  in2.  The  total 
exposed  area  (A)  =  44.66  +  28.274  +  2.828  +  9.425  «  85.19  in2. 
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Fig.  F-4  Nozzle  surface  area  as  a  function  of  nozzle  throat  diameter  for 
BATES  motor  nozzles . 


